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A Convex Sum-of-Squares Approach to Analysis, 
State Feedback and Output Feedback 
Control of Parabolic PDFs 

Aditya Gahlawat and Matthew. M. Peet 


Abstract —We present an optimization-based framework for 
analysis and control of linear parabolic Partial Differential 
Equations (PDEs) with spatially varying coefficients without 
discretization or numerical approximation. Eor controller syn¬ 
thesis, we consider both full-state feedback and point observation 
(output feedback). The input occurs at the boundary (point 
actuation). We use positive definite matrices to parameterize 
positive Lyapunov functions and polynomials to parameterize 
controller and observer gains. We use duality and an invertible 
state variable transformation to convexity the controller synthesis 
problem. Einally, we combine our synthesis condition with the 
Luenberger observer framework to express the output feedback 
controller synthesis problem as a set of LMI/SDP constraints. 
We perform an extensive set of numerical experiments to 
demonstrate accuracy of the conditions and to prove necessity 
of the Lyapunov structures chosen. We provide numerical and 
analytical comparisons with alternative approaches to control 
including Sturm Liouviile theory and backstepping. Einally we 
use numerical tests to show that the method retains its accuracy 
for alternative boundary conditions. 

Index Terms —Distributed parameter systems, partial differen¬ 
tial equations (PDEs), control design, sum of squares. 


I. Introduction 

Partial Differential Equations (PDEs) are used to model 
quantities which vary in both space and time with early 
examples including the D’Alembert wave equation (1746); 
the Euler-Bernoulli beam (1750); the Euler equations (1757); 
and the Eourier heat equation (1822). Today, the use of PDE 
models has expanded to include phenomena such as the mag¬ 
netohydrodynamics of plasma in a fusion reactor ||43]| . tumour 
growth, infectious diseases, and ecological succession ll25l 
Chapter 11]. However, despite the variety of phenomena 
modeled by PDEs, compared to the literature on Ordinary 
Differential Equations (ODEs), our knowledge of how to 
analyze and control PDEs remains incomplete. 

Consider the following class of scalar-valued anisotropic 
parabolic PDEs with input u{t) G R, 

Wt{x,t) = a{x)wx^{x,t) -I- b{x)w^{x,t) + c{x)w{x,t), (1) 
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X G [0,1], t > 0, which has output v{t) = w{l,t) G R. and 
mixed boundary conditions of the form 

w{0,t)=0, Wxil,t) = u{t), (2) 

where a, b and c are polynomials with a{x) > a > 0, for 
X G [0,1]. We assume the controller is parameterized by scalar 
Ri and function R 2 as u{t) = Riw{l,t)+ R 2 (x)w(x,t)dx 
where w is an estimate of w obtained from some set of 
observer dynamics. The objective of the paper is to propose 
an optimization-based method for determining controller gains 
Ri and i ?2 and observer dynamics which minimize certain 
closed-loop gains. 

Control of PDE models is a challenging problem in that 
slight variations in the type of PDE, boundary conditions, etc. 
may dramatically alter properties of the solution E 2 I . The 
model defined above is classified as an anisotropic parabolic 
PDE with point inputs and point outputs. The term anisotropic 
means that the values of the coefficients a{x),b{x) and 
c{x) depend on the spatial variable x G [0,1]. Examples 
of anisotropic systems include heat conduction with non- 
homogeneous conductive properties or a wave propagating 
through a medium of varying density. The term point input 
(boundary actuated) means that the control input determines 
one of the boundary values and therefore has no direct mea¬ 
surable effect on Equation ([T]i. This is in contrast to the case 
of distributed inputs, wherein the control effort is spread over 
some measurable subset of the domain. In a similar manner, 
the term point output means that the sensor measures the state 
at a single point in the domain and hence the output operator 
is unbounded in the L 2 induced norm. 

Perhaps the most common approach to analysis and control 
of PDEs is based on the use of discrete approximation. Such 
approximation techniques typically use a model reduction 
wherein the PDE is approximated by a set of ODEs. Einite- 
dimensional linear control theory is then used to analyze 
stability and design control laws for the finite-dimensional 
approximations ll2^ . ll2Tl . Eurthermore, results have been 
obtained which show that as the order of the discrete approx¬ 
imation increases, stability of the closed-loop approximations 
will eventually imply stability of the closed-loop PDE. A 
disadvantage of the discrete approach, however, is that the 
required order of the approximation cannot be established a 
priori. Consequently, the stability of any particular approx¬ 
imation is not guaranteed to imply stability of the actual 
PDE. Eor this reason, among others, there has been some 
interest in finding approaches to analysis and control which 
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can be applied directly to the PDE model without the use of 
discretization or numerical approximation. Such methods are 
sometimes termed direct or infinite-dimensional. 

There has been signihcant progress in the use of direct 
methods for control of PDE systems. One approach is to 
express the control problem as the solution to a set of 
operator-valued Riccati equations. This approach was applied 
to distributed input/distributed output optimal control problems 
in ED . The problem of point actuation with full-state feedback 
was considered in II 2 TII (and related work) and extended in ll20l 
to output feedback controller synthesis through the use of a 
Luenberger observer. An alternative Riccati-based approach 
for static output feedback of a certain class of well-posed 
operators can be found in ll35l . 1^ . ll42l . A limitation of these 
Riccati-based methods, however, is that they rely on finite¬ 
dimensional numerical methods for obtaining the operator¬ 
valued solution. While convergence of these approximations 
has been demonstrated ED, for a given level of approxima¬ 
tion, it is not possible to determine whether existence of a 
solution implies the closed loop is stable when applied to the 
original PDE. 

Backstepping lITSl is a popular and well-developed method 
for boundary control of parabolic PDE systems. This approach 
is based on the use of a boundary controller to transform 
the PDE to a simpler model for which the existence of a 
decreasing Lyapunov function has previously been established. 
The backstepping approach is commonly used in the literature 
and has been extended to many classes of PDE systems - see, 
e-g- mi, EH, im. Ell. A highlight of the backstepping 
method is that for certain types of system, stabilizability 
guarantees the existence of a backstepping transformation. 
However, a drawback of the backstepping approach is that 
it is not based on optimization, but rather typically requires 
numerical integration of a PDE in order to obtain the stabi¬ 
lizing controller - thereby making extensions to robust and 
optimal control more difficult. Although a complete survey 
the of the literature on direct control of PDEs is beyond 
the scope of this paper, we do note some other significant 
results on the use of Lyapunov functions for analysis and 
control of infinite dimensional systems including: a rotating 
beam 0; quasilinear hyperbolic systems 0; and control of 
systems governed by conservation laws 0. As an alternative 
to Lyapunov-based methods, a classical spectral approach to 
stability and stabilization is based on Sturm-Liouville theory. 
In particular, the differential operators which dehne the PDEs 
in this paper can be adapted to the Sturm-Liouville framework, 
from whence one can attempt to determine stability and 
design static output-feedback controllers. As is demonstrated 
in Section IXIIII however, the use of dynamic output feedback 
offers considerable advantages over this classical framework. 

The goal of this paper is to design stabilizing static state 
feedback and dynamic output feedback controllers for PDE 
systems. Our approach is inspired by the use of Linear Matrix 
Inequalities (LMIs) and Semi-Definite Programming (SDP) 
in control of ODEs. Eor stability analysis, as discussed in 
Sections |V]and|Vll we use positive definite matrices to create a 
linear parametrization of a cone of Lyapunov functions which 
are positive on the Hilbert space L 2 - Specifically, the Lyapunov 


functions have the quadratic form V = {Z(w), PZ(w)} 
where w G L2 is the infinite-dimensional state, P is a positive 
definite matrix and Z is a fixed vector of multiplication and 
integral operators with monomial multipliers and kernels. The 
derivative of the Lyapunov function is likewise constrained 
to be a negative definite quadratic form. If such a Lyapunov 
function exists it directly proves stability of the PDE - i.e. there 
is no numerical approximation. Eor state-feedback controller 
synthesis, the controller, as defined above, is parameterized by 
a scalar Ri and a function i? 2 - Combining these gains with the 
quadratic Lyapunov functions used for stability analysis yields 
synthesis conditions which are bilinear in the design variables. 
However, as described in Sections fVIIl and IVIIII by dehning an 
invertible state transformation and a variable substitution, we 
derive synthesis conditions which are linear in the optimization 
variables. Next, in SectionH^we introduce a class of infinite¬ 
dimensional Luenberger observers with observer gains, again 
parameterized by the coefficients of polynomials. Again, using 
the Lyapunov function from Section |V] and the invertible state 
variable transformation from Section IVIII we obtain SDP- 
based observer synthesis conditions. Einally, in Section]^ we 
verify the accuracy of the method with a series of numerical 
tests which indicate that the proposed stability conditions are 
accurate to several decimal places and suggest that for any 
suitably controllable and observable system, the algorithm 
will return an observer-based controller. This is followed by 
Section IXIIII wherein we include numerical and analytical 
comparisons with other results in the literature, including 
Sturm-Liouville and backstepping. 

A signihcant contribution of the paper, in addition to a 
new approach to analysis and control of PDEs, lies in the 
Hexibility of the optimization-based approach. Specihcally, as 
the use of LMIs for control of ODEs enabled the held of 
robust control, so too does our LMI/Lyapunov-based approach 
to control of PDEs allow the extension to analysis and control 
of PDEs with parametric uncertainty, PDEs with nonlinearity, 
multivariate PDEs and PDEs coupled with ODEs or delays. 
Einally, we note that our approach is complementary to several 
recent results in the use of LMIs for stability and control of 
PDEs, including, e.g. our early work in ED, modeling and 
control of nonlinear dynamic systems in E9l , stability analysis 
of semilinear parabolic and hyperbolic systems in ifT^ and the 
numerous results contained in ll26l . 

II. Notation 

We denote the vector space of m-hy-n real matrices by 
R™x" and the subspace of symmetric matrices by §" C 
Yvjjgj-g the multiplicative and additive identities are 
denoted by G S" and Om,n G respectively. Eor 

P G S", P > 0 (P > 0) denotes that P is a positive 
dehnite (positive semi-dehnite) matrix. The spaces of n—times 
continuously differentiable and inhnitely differentiable func¬ 
tions on an interval PE C K are denoted by C'^{W) and 
C^iW), respectively. In a similar manner, C'"’™(Wi,W 2 ) 
represents the space of n and m—times continuously dif¬ 
ferentiable functions on intervals Wi C R and W2 C R, 
respectively. The shorthand Ux and ut denote the partial 
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derivative of u with respect to independent variables x and 
t, respectively. For a bivariate function, f{x,y), we denote 
Dif := fx and £> 2 / := fy - i.e. Di is differentiation with 
respect to the first variable and D 2 is differentiation with 
respect to the second. In a similar manner, := f^x and 
D\ := jyy. Recall L 2 {W) is the standard Hilbert space 
of square Lebesgue integrable functions with standard norm 
and inner product. We use H^{W) to denote the Sobolev 
subspace H'^{W) := € L 2 {W) : ^ S £ 2 (W^)| with 

inner product (x,t/)^„ = El=o We occa¬ 

sionally let £ 2 ( 0 ,1) := £ 2 ([ 0 , 1]) and i£"(0,1) := £f”([0,1]). 
For normed spaces X and Y, C{X, Y) denotes the Ba¬ 
nach space of bounded linear operators from X to y with 
induced norm ||C/||£ := sup|| 2 ,||^^ 2 ll^^a;||F and we denote 
C{X) ■= C{X,X). We define Zd{x) £ to be the 

column vector of all monomials in variables x of degree d 
or less arranged in increasing lexicographical order. We often 
use the notation Zd{x,y) := Zd{[x\y\) to denote the vector 
of monomials in both x and y. For any function £ G £2 
we use A4t ■ L 2 ^ L 2 to denote the multiplier operator 
defined by T. i.e. {A4tw){x) = T{x)w{x). For any functions 
M,Ki,K 2 £ we define 

= M{x)w{x)+ f Ki{x,^)w{()d^+ [ K 2 ix,^)w(()d^. 

JQ Jx 

(3) 

III. Problem Statement 

For the system of Equations ([T]i - dU, the strict positivity of 
a{x) implies that the differential operator defining the PDE 
is uniformly elliptic uni Section 6.1]. This means that w 
diffuses from higher density to lower density, a property which 
is representative of most physical systems. The choice of 
sensor and actuator location is somewhat arbitrary. Eor the 
heat equation, input Wx{i,t) = u{t) would represent heat 
flow into the rod and the output v{t) = tii(l, t) represents the 
temperature of the rod at that point. Note that the results of 
this paper can be adapted to Dirichlet, Neuman and Robin 
boundary conditions with only slight modifications to the 
conditions and proofs. These extensions are addressed in 
Section IXlV] 

The goal of this article is to design algorithms which resolve 
the following problems: 

1) Stability Analysis: Establish global exponential stability 
of the trivial solution w = 0 of the autonomous system 
u{t) = 0 and determine the exponential rate of decay S. 

2) State feedback control: If the autonomous system is 
unstable, construct gains i?i G K and R 2 {x) £ C°^(0,1) 
such that if 

u(t) = Iliw(l,t) + f R 2 {x)w{x,t)dx, (4) 
Jo 

then the trivial solution m = 0 is globally exponentially 
stable with some desired rate of decay, y. 

3) Output feedback control: If only output feedback is 

available {v{f) = ti;(l, t)), construct gains £1 G 1) 


and £2 G K. such that for stabilizing gains i?i and R 2 , if 

u{f) = Riwfi^f) + / R 2 {x)w{x,t)dx, 

Jo 

where w satisfies 

Wt(x, t) =a{x)wxx{x,f) + h{x)wx{,x, t) 

+ c{x)w{x, t) -I- £i(a;) {v(t) - v{t)), 

(5) 

for v{t) = w{l,t) and v{t) = w{l,t) with boundary 
conditions 

u;(0, t) = 0, Wxil, t) = u(t) + £2 (w(f) - v(t)), 

(6) 

then the trivial solution m = 0 of Equations ([T]) - (|2]l is 
globally exponentially stable. 

Note that if we consider only bounded linear operators, then 
the structure of the controller in Q is not restrictive, as any 
bounded linear functional can be represented in this way using 
only the integral form (second term). However, we also would 
like to consider unbounded operators and hence we include 
the term Riw{l,t) as well. If controllers of this form prove 
inadequate, then one can generalize the structure further to 
include terms such as R3(x)wx(x,t)dx as in IfTSl . 

The choice for the strucmre of the Luenberger observer was 
similarly determined in an ad-hoc manner through inclusion of 
terms necessary to achieve separation of controller synthesis 
and observer design objectives. That is, the goal of the 
observer is to stabilize the dynamics of the estimation error 
e = w — w and the terms in Equations (ID - (|6]l were chosen 
as the minimal necessary to achieve this objective. Again, 
this structure mirrors the structure of observers found in the 
backstepping approach. 

A. Existence and Uniqueness 

We now briefly discuss the uniqueness and existence of 
solutions. Define the operator 

(f d 

A = a{x)-^ + b{x)— + c{x). (7) 

dx^ dx 

It is known that the operator A restricted to space 

Vo = {w £ H\0A) : u;(0) =«;,(!) =0}, (8) 

generates a strongly-continuous semigroup, or a Cg- 
semigroup, on £ 2 ( 0 , 1 ) (see, e.g., ||7] Section 2.1]). More 
precisely, one can represent A as the negative of a Sturm- 
Liouville operator on Vq and hence, using the spectral prop¬ 
erties of a Sturm-Liouville operator, it can be proven that A 
restricted to Vq generates a Cg-semigroup on £ 2 ( 0 , 1 ) Il9l . 
Thus, using Theorems 3.1.3 and 3.1.7 in 121 we conclude that 
in the autonomous case (u(<) = 0), for any initial condition 
wq £ Vq there exists a unique classical solution of dTJ - (ID. 

Eor the state-feedback case, using a fixed point argument 
similar to the one presented in ID it can be shown that for 
i?i G K and R 2 £ £ 00 (0,1), the closed loop system dD - (O 
with ^ 

u(t) = Riw(l,t) + / R2(x)w(x, t)dx, 

Jo 
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admits a unique local in time solution w G T), [0,1]), 

for T > 0 sufficiently small, for any initial condition wq G V, 
where 

V = {wG 1): w(0) = 0 and 

r(;a^(l) = + / R2{x)w{x)dx}. 

Jo 

(9) 

Thus if we can establish that any solution of the closed 
loop system decays exponentially, then this implies the ex¬ 
istence and uniqueness of a unique classical solution w G 
oo), [0,1]) for any wq G T>. The proof of this state¬ 
ment has been omitted, but follows the arguments presented 
in ||2] Section 6]. 

Finally, consider the observer-based controller as defined in 
Equations ([T]i - (|2|i and 0 - (|6l). Define the estimator error as 
e = w — w, which is governed by 

etix,t) = a{x)exx{x,t)+b{x)exix,t)+c{x)e{x,t)+Li{x)e{l, 

( 10 ) 

with boundary conditions 

e(0,f) = 0, e^{l,t) = L 2 e{l,t). (11) 

It has been established in ifT^ Section 2] that for Li G 
(^^(OjOo) and L 2 G K., Equations (fTOl i - (fTTT i. if expo¬ 
nentially stable, admit a unique local in time solution e G 
T), [0,1]), for T > 0 sufficiently small, for any initial 
condition cq G De, where 

Ve = {w € iF^(0,1): tt;(0) = 0 and iCa;)!) = ^ 21 ^( 1 )}. 

( 12 ) 

Therefore, if we can establish that any solution of the 
coupled closed-loop dynamics decays exponentially, then the 
local in time solution can be extended to a classical solution 
e G 00 ), [0,1]) for any initial condition Cq G 'De- 

IV. A Eramework for Stability Analysis and 
Control 

Our approach is motivated by the use of LMIs for optimal 
control of finite-dimensional systems. Eor example, consider 
the autonomous finite-dimensional ODE 

x{t) = Ax(t), 

where x{t) G K.". This ODE is exponentially stable if and 
only if there exists a positive definite matrix P G §" such that 

+ PA<0. 

Eeasibility of this LMI implies that the Lyapunov function 
V{x) = x"^Px is positive definite and its derivative along 
solutions V{x) = x'^{A^P + PA)x is negative definite. 
Eor stability of PDEs, our approach is to use positive ma¬ 
trices to define positive quadratic Lyapunov functions, ex¬ 
cept that instead of V{x) = x^Px, we will use the form 
V{w) = {Z(w), PZ(w)), where Z : L 2 —>■ is a vector 

of bases for a subspace of linear operators on L 2 (similar 
to how X = [xi,--- is a vector of bases for the 

space of linear functions on R"). In our case, however, Z 
parameterizes a subspace of multiplier and integral operators 


with polynomial multipliers and semi-separable kernels. Then, 
if P > el, it has a symmetric square root and hence 
V{w) = (^P^Z{w), P^Z{w)'^ > e||w|p. Eor the time 

derivative, we will similarly require V{w{t)) + fiV{w{t)) = 

— {Z{w{t)),QZ{w{t))), for some scalar p > 0 and Q > 0 

where here and throughout the paper we denote by V the 
function which satisfies L(ui(f)) := for any 

solution of the associated PDE - i.e. the derivative along 
solutions or time-derivative. Existence of such P,Q>0 
implies exponential stability of the system. As was done for 
LMIs in finite-dimensional systems, this approach can then 
be extended to controller and observer synthesis, as outlined 
below. 

a) Controller Synthesis: Eor controller synthesis, again 
consider the LMI approach for the finite-dimensional system; 

x{t) = Ax{t) + Bu(t), 

where x{t) G MA and u(t) G R"*. Eor this system, there 
exists a stabilizing state feedback controller of the form u{t) = 
Rx{t) if and only if there exists a positive definite matrix P 
and Y G such that 

{AP + BY) + [AP + BY)^ < 0. 

If this LMI is feasible, then for R = YP~^, the Lyapunov 
function V{x) = x'^P~^x is positive definite and has time 
derivative 

V = X^(P-^A + P-^BR + A^P-^ + (PP)^p-i)x 
= {P-'^X)'^{AP + BRP + PA^ + {BRP)'^){p-^x) 

= y^iAP + BY + PA^ + {BY)'^)y < 0, 

where y = P' ^X. The extension of this LMI approach to 
PDEs is to search for a positive definite operator P = Z*PZ 
for some P > 0 and operator y, defined by {yz){z) := 
yiz(l) + fg Y 2 (x)z(x)dx, such that if m = TZw = yV~^w, 

u{t) = Riw{l,t) + / R 2 {x)w{x,t)dx 

Jo 

= Yi{V~^w){l,t)+ ( Y 2 {x){P~^w){x,t)dx, 

Jo 

then the Lyapunov function V = (w,V~^w) = 

(^Z{P~^w),PZ{P~^w)) satisfies V{w{t)) + 2yV{w{t)) = 

— (^Z (V~^w) ,QZ (V~^w)) for some scalar p, > 0 and 
Q > 0, which implies the closed-loop system is exponentially 
stable. This is detailed in Section IVIIII 

b) Observer Synthesis: As mentioned previously, for 
observer design, we use a Luenberger observer and a sep¬ 
aration principle to decouple the error dynamics as defined 
in Equations (fTol) - (fTTI) . Eor a finite-dimensional Luenberger 
observer, where the output is v{t) = Cx{t), the estimator 
dynamics are defined using the controller gain F and observer 
gain L as 

X = {A + LC)x — Lv(t) + Bu{t). 

If u{t) = Fx(t), then the error dynamics become 
e(t) = {A + LC)e{t). 
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Existence of an observer gain which renders the error dynam¬ 
ics stable is equivalent to the existence of a P > 0 and T such 
that 

PA + TC + A^P + C^T'^ < 0. 

If this LMI is feasible, then for L = P~^T, the Lyapunov 
function V (e) = e^Pe is positive definite and has derivative 

V(e) = e^{PA + PLC + A^ P + C^L^P)e 
= e^{PA + TC + A^P + C'^T'^)e < 0. 

For the infinite-dimensional PDE, we have two observer gains 
which we construct as 

Liix)=p-\mx)+T3ix)) and L 2 =p-\T 2 ), 

for some gains Ti, T 2 and T 3 and where V = Z*PZ for 
some P > 0. We then use the Lyapunov function V{e) = 
{Z{e)^PZ{e)) and search for a Q > 0 such that E(e) = 
— {Z{e),QZ{e)) < —SV(e) < 0, for some 5 > 0. This is 
detailed in Section |I2 


V. Sum-of-Squares Lyapunov Functions with 
Semi-Separable Kernels 

In this Section, we define the map Z and show how this map 
is used to construct Lyapunov functions of the form V{w) = 
{Z(w), PZ(w)). This approach is based on prior work, as 
described in ||28l. Specifically, we define 


{Zw)ix) = 


Zd^{x)w{x) 

J^Zd^{xA)w{0d^ 
Jo ^d 2 (x,C)w( 0 di 


where recall Zd^ (x) and Zd 2 (x, are the vectors of all 
monomials of degree di and d 2 or less, starting with 1 . 

Theorem 1. Given di,d 2 G N and e > 0, e £ M, let 
Zi{x) = Zdi(x) and Z 2 {xA) = Zd 2 (xA), with n = di + 1 
and m = ^{d 2 + 2 ){d 2 +l) denoting the length of these vectors, 
respectively. Suppose that there exists a matrix P £ 
such that 


P = 


Pii - 


e 

n — 1,1 On — l,n — 1 

P12 

Pl 3 

X12 

P22 

CO 

Xl 3 

-^23 

P33 


> 0 , 


(13) 

where Pij is a partition of P such that Pn £ §"',P 22 £ 
and P 33 £ Now let 

Mix) = Zi(x)^PiiZi(cr), (14) 

KiixA) = Ziix)^P,2Z2ix,0 + Z2i^,x)'^P3iZ,i0 

+ [ Z2iri,x)'^PooZ2{r]A)d'q 
+ [ Z2ir],x)'^P32Z2ir],0dr] 


+ / Z2ir],x)'^ P22Z2ir]A)dv, 

J X 

K 2 ix ,0 = Ki{^,x). 


(15) 

(16) 


Then 


Viw) = / wix)Mix)wix)dx 
Jo 

+ f f w{x)Ki{xA)wi^)d^dx 

Jo Jo 

+ f ( w{x)K2ixA)'wi^)d^dx (17) 
Jo Jx 

= {Z{w),PZiw)) = (^P^Ziw), P^Ziw)^ > e||u;|p. 


Proof: The proof follows directly from the definition of 
Z and the Sum-of-Squares representation of V. ■ 

The form of the Lyapunov function defined by Theorem[T]in 
Equation (fTTI i is somewhat atypical for the study of parabolic 
PDFs. A more commonly used version would be V(w) = 
w(x)M(x)w(x)dx or even yet V(w) = w(x)Mw(x)dx 
for M > 0. Such forms can be obtained as a special case 
of Theorem [T] when Pij — 0 for i j ^ 1. However, as 
we discuss in Section IXIIl neglect of the Ki and K 2 terms 
results in significantly less accurate conditions for stability and 
control. 

For polynomials M, Ki and K 2 , let X{m,Kx,K 2 } be de¬ 
fined as in Q. If M, Ki and K 2 satisfy the conditions 
of TheoremlH then Viw) = {w, X^m,Ki,K 2 }'>x) > e||tp|P, 
which implies the operator positive definite 

and furthermore, coercive. Moreover, since M, Ki and K 2 are 
polynomials, the operator is bounded, which implies that there 
exists a 6 > 0 such that e||w|P < Viw) < 0||w|p. Finally, 
the constraint (fTbl l in Theorem [T] implies that the operator 
X{M.K^,K 2 } is self-adjoint. 

As discussed in Section HYl Theorem [T| allows us to use 
positive matrices to parameterize positive Lyapunov functions 
of the Form (fTTl i. By expanding these forms, the coefficients 
of the polynomials M, Ki and K 2 are linear combinations 
of the elements of P > 0. Furthermore, if we can express 
the derivative V in the Form ( fTTI i, where the coefficients are 
again linear combinations of the elements of P, then we can 
enforce negativity of the derivative along the solutions w by 
using Viw) = — {Ziw),QZiw)) to equate these coefficients 
to those defined by Q > 0. Constructing the matrices which 
relate the elements of P and Q can be automated using 
MATLAB toolboxes for polynomial manipulation such as 
MULTIPOLY, contained in the package SOSTOOLS ll29l and 
further developed in our package DELAYTOOLS Il28]l . 

For polynomials M, Ki and K 2 , we represent the constraint 
{w,X{M,Ki,K 2 }'dt) = (Ziw), PZiw)) for some P > 0 as 
{M,Ki,K 2 } £ '^{d^,d 2 ,e} where 

S{di,d 2 .e} ■={M,Ki,K 2 : M,Ki,K 2 satisfy 

Theorem [T] for di, d 2 , e}. 

The constraint {M, Ki, K 2 } £ Zj^di,d 2 ,e} ™ LMI constraint 
in the coefficients of the polynomials M, Ki and K 2 and 
the unknown matrix P > 0. In this way, the shorthand 
{M,Ki,K2} £ 'E.[di,d2,e} allows us to define LMI constraints 
implicitly. 
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VI. A Test for Stability 


the classical solution w of (HI - (0 satisfies 


In this section, we use the results of the previous section 
to test the existence of a Lyapunov function which establishes 
stability of the scalar parabolic PDE defined in Equations ([T]i - 
dill. Recall the autonomous {u{t) = 0) form of the PDE 

wt{x,f) = a{x)wxx{x,t) + b{x)wx{xfi) + c{x)w{x,t), 

( 18 ) 

■u;(0,f)=0, u; 3 ;(l,f) = 0. (19) 


The main technical contribution of this section is reformulating 
the derivative of the Lyapunov function V in (fTTI) in the form 
of Equation (fTTl i. This is achieved in the following theorem 
wherein we obtain functions M, Ki and K 2 such that 


V{w) < / w{x)M{x)w{x)dx 
Jo 

pi pX 


0 ^0 


w{x)Ki {x, ^)w{^)d^dx 

ft /•! 


w{x)K 2 {x, ^)w{Cjd^dx. 


0 J X 


Note that the inequality in this expression is deliberate, i.e., 
certain negative semidehnite terms have been left out of M, 
Ki and K 2 . 

Before giving the main theorem, we define the following 
linear map, fig, which relates functions M, Ki and K 2 to an 
upper bound on the time-derivative of the Lyapunov function 
defined by these functions. Specifically, we say that 

{M,K,,K2}~^s{M,K„K2), (20) 


if 




—a{x)M{x) — b{x)M{x) 


+ 2 —[a{x){Ki{x,^)-K 2 ix,^))] 

TT^ 

+ 2M {x)c{x) — 


J ^=x 

( 21 ) 


Kiix ,0 = 


dx 


dx 


[a{x)Ki(x,0] - Hx)Ki{x,C) 


+ 


d 


d 


[a{OKi{xM-b{OKi{x,0 


+ {c{x) + c{^)) Ki{x,^), (22) 

K2ix,0=Kii^,x). (23) 


Theorem 2. Suppose that there exist scalars e,d > 0, 
di,d 2 -,di,d 2 G N and polynomials M, Ki and K 2 such that 

{M,Ki,K2} G Sdi.da.ei 

{_M - 25M, -ki - 2SKi, -k2 - 2 SK 2 } G 

( 6 ( 1 ) - a,(l))iTi(l,a:) - a(l)(f9iiTi)(l, x) = 0 , 

( 6 ( 1 ) - a,(l))M(l) - a(l)M,(l) < 0 , 

K2{0,x) = 0, 

where {M,Ki,K 2 } '■= flsiM, Ki, K 2 ). Then for any initial 
condition r(;(0) G Do> there exists a scalar 7 > 0 such that 


lk(f)|| < 7lk(0)||e t > 0, 

where Dq is defined in Equation (©. 

Proof: Recall the operator defined 

in (O. As discussed in Section |III] for any w{0) G Do the 
autonomous system admits a unique classical solution. By 
Theorem [U if {M, Ki, K 2 } G ^ 2,0 then 

V{w) = {w,X[m,Ki,K 2 }U>) = {w.Vw ), 

satisfies e||tu||^ < V{w) < for some 0 > 0. The 

calculation of the time derivative V and its reformulation is 
lengthy. It involves integration by parts, the Wirtinger inequal¬ 
ity and the assumption a(x) > a. Eor this reason, we have 
included this proof in the appendix as Lemma [3 Continuing, 
by Lemma |3 for any w which satisfies Equations (fTSl l - ( fT^ , 

V(w{f)) < (w(i)>‘^{M,iVi.iV 2 }tt'( 0 ) • 

Now, since {-M - 26M, -Ki - 2 SKi,-k 2 - 28 X 2 } G 
=dj.J 2 , 0 ’ we have that ^ and thus 

— V{w) — 26V{w) > 0. This implies that ■^V{w{t)) < 
—2SV{'w{t)) for all f > 0. Thus, V{w{t)) < L(u;(0))e“^'^*. 
Concluding, we have that 

||u;(f)|| < 7 ||w( 0 )||e“‘^‘, 7 = 

■ 

Note that using the arguments in the proof of Q Theo¬ 
rem 5.1.3], the above result holds for weak/mild solutions 
where the initial condition need only satisfy wq G ^ 2 ( 0 ,1). 

To test the conditions of Theorem |2] the variables are 
the coefficients of the polynomials M, Ki and K 2 . The 
coefficients of M, ki and k 2 are then linear combinations 
of these variables. Einally, the constraints G d 2 .e ^6 
LMI constraints, as discussed in Section |V] Constructing the 
matrices which map these coefficients can be automated using 
SOSTOOLS or DelayTOOLs. The algorithm used can be 
adapted from the algorithm presented for output feedback 
controller in Section |X] Application of the conditions of 
Theorem |2] to several numerical examples can be found in 
Section 

VIE Inversion and State Transformation 

As discussed in Section |IV] for controller synthesis, we 
will use a state variable transformation z = V~^w so that 
{Zh, PZ{V~^w)) = (6,, w). Define V = X{m,Ki,K 2 }^ where 
^{m,Ki,K 2 } defined in (|3ll. Then V has the form 

{Pz){x)^{x)z{x)+ f Ki{x,^)zk)d^+ 

Jo 

where if {M,Ki,K 2 } G the operator is coercive 

with {w,Vw) > e||tx|p. Operators of this type are a com¬ 
bination of a multiplier operator and two integral operators. 
Eurthermore, since Ki and K 2 are polynomials, there exist 
polynomials E) and Gi such that Ki{x,^) = Fi{xY'Giff) 
and K 2 {x,^) = E 2 (x)^G 2 (^). This implies that the two 


/ K 2 ix, 0 zi 0 d^, 
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integral operators can be combined into a single integral 
of the form K(x,^)z(^)d^ where K is a kernel of the 
semiseparable type. That is, there exist functions Fi and Gi 
such that 

’ otherwise 

Integral operators with semiseparable kernels are used to repre¬ 
sent the input-output map of well-posed Linear Time-Varying 
(LTV) systems, as explored in IfT^ Section 1.4, Theorem 4.1]. 
These operators have certain properties which make them well- 
suited for use in Lyapunov functions. Specifically, they are not 
trace-class, which means that their eigenvalues may not be 
summable. Moreover, as discussed in ifTbl Section 11.2], since 
M{x) > e > 0, 'P~^ is a bounded linear operator and can 
be calculated explicitly, as in the following theorem, which is 
adapted from ifT^ Section 11.3, Theorem 3.1]. 


Theorem 3. Suppose that {M,Ki,K 2 } G Syi.da.e} far 
some di,d 2 ,£ > 0 with Ki{x,^) = F{xyG{^) and 

K2 {x,C) = G{x)'^F{faj. Let V G £(^ 2 ( 0 ,!)) be defined as 
'P — where X{m,Kx,K' 2 } is as defined in Q. 

Define 


B{x) 


'G{x) 
F{x)\ ’ 


G{x) = [F{x)^ 


H =[Ni+N2U{l)]~^ 


Ni 


I 0 
0 0 ’ 


N2 = 


N2U{1) 

0 o' 

0 I ’ 


-GixV] , 


and U{x) = lim„^oo Un(x), where 

t7„+i(x) = 1- r B{OM{(,r"G{OUn{m. (24) 
fa 

and Ui = L Then, the inverse of the operator V is given by 


(V ^w) {x) =M{x)w{x) + / K^{xfa)w{fa)d^ 

Jo 

+ [ K2{x,0wi0d^, 

J X 

M{x) =M{x)~^, 


K,{x,0 =M{x)-^C{x)U{x){H - I)U{0~"B{0Mi0-^ 
K^ix,0 =M{x)-^C{x)Uix)HU{0~"Bi0Mi0~". 


Note that since M(^) > e, is bounded and 

continuous and hence the matrix of rational functions 
B{^)M{^)~^G{^) is bounded and continuous. Therefore, it 
follows from E Chapter 3] that the uniform limit U{x) 
exists and is non-singular for x G [0,1]. Since U{x) is non¬ 
singular on [0,1], the matrix H is well defined. Therefore, by 
construction M,K]^,K 2 G G°°. Furthermore, note that since 
V satisfies e||w||^ < {w^Vw) < for some 0 > 0, then 

l/ 0 ||ru|p < (w,V~^w) < l/e||r(;|p. 

Theorem |3] not only proves existence, but gives a practi¬ 
cal method for constructing the state variable transformation 
V~^ for which (^Zh, PZ('P~^w)^ = {h,w). Specifically, if 
we truncate the sequence Un and approximate M(x)~^ by 
a Chebyshev series, then construction of the functions M, 
Kj and K 2 is simply a matter of polynomial multiplication 



and integration, which can be performed in MATLAB or 
Mathematica. In practice, we have found that [/„ converges 
after only a few iterations. To illustrate, in Figure [T] we have 
applied this approach to a given {M,Ki,K 2 } G Si 11 and 
plot ||w — Wf^pwW as a function of n for the arbitrarily 
chosen function w{x) = x(x—QA){x—l). Here denotes 
the construction for defined in Theorem [3 with U{x) 
replaced by Un+i{x). In this case, n = 5 yields an L 2 norm 
etTor of ~ 10“®. In this example, we approximated 
using the first five terms of its Chebyshev series. 

Finally, we emphasize that construction of 'P~^ is not part 
of the optimization algorithm, but rather is performed after the 
algorithm has solved the controller synthesis problem (to be 
defined in the following section) and returned the polynomial 
variables M, Ki, and K 2 . 

VIII. State-Feedback Controller Synthesis 
Our approach to controller synthesis is based on the use 
of a state variable transformation y — V~^w which, by 
Theorem [3l is guaranteed to exist for any V = X^m,Ki,K 2 } 
defined by {M,Ki,K 2 } G E[di,d 2 ,e}- Specifically, we will 
use the Lyapunov function V{w) = (V~^w,w') = {y,Vy). 
Ignoring the input for the moment and using the operator A 
defined in Equation (|2l), the time-derivative of this function 
yields the dual stability condition 

V[w) = 2 {V-^w, Aw) = 2 {y, AVy) < 0, 

which we must enforce for all j/ G £ 2 - The critical point 
is that the operator 'P~^ does not appear explicitly in the 
stability condition. Rather its existence is only inferred from 
the constraint on V that {M,Ki,K 2 ] G The next 

step in our approach is to combine this dual stability condition 
with a variable substitution through the use of a controller of 
the form 

fat) = Yi{V~^w){l,f) + f Y 2 {x){'P~^w){x,t)dx 
Jo 

= Riw{l,t) + / R 2 {x)w{x,t)dx, 

Jo 

wherein we have replaced the original controller gains Ri and 
i ?2 with the new variables Yi and 12 - Once Yi and Y 2 are 









determined by the SOS solver, the actual gains i?i and R 2 
can be recovered by computing 'P~^ and applying the formula 
listed here. 

Before giving the main theorem, we recall that the input 
enters the dynamics as 


wt{x,t) = a{x)wxx{x,t) + b{x)wx{x,t) + c{x)w{x,t), 


w{ 0 ,t)= 0 , Wxil,t) = u{t). 


(25) 

(26) 


The goal, then, is to dehne conditions on P (which dehnes 
M, Ki and K 2 ) as well as on Yi and the polynomial Y 2 such 
that the closed-loop system is exponentially stable. 

To simplify exposition, we now define the following linear 
map, Oc, which relates functions M, Ki and K 2 to an upper 
bound on the time-derivative of the Lyapunov function dehned 
by these functions for the controller dynamics. Specihcally, we 
say that 

{M,K,,K 2 } ■.= n,iM,K,,K 2 ), 


(27) 


if 


M{x) = {axx{x) - bx{x))M{x) + b{x)Mx{x) 

TT^ 

+ a{x)Mxx{x) + 2c{x)M{x) - 

+ a{x) 2-^[Ki{x,C) - K 2 {x, 0 ] 

ki{x,C) = a{x){DlKi){x,C) + Kx){DiKi){x,C) 


(28) 


amDlK^){x ,0 + b{ 0 {D 2 Ki){x ,0 

+ {c{x)+c{^))Ki{x,^), (29) 


k2{x,0 = ki{^,x). 


(30) 


Theorem 4. Suppose that there exist scalars e, p, > 0, 
di,d 2 ,di^d 2 G N and polynomials M, Ki and K 2 such that 

{M,Ki,K 2 ) and K 2 {Q,x)=Q. 

Further suppose 

{-M - 2/rM, -K, - 2pKi, -K 2 - 2 ^X 2 } S 

where {M,iTi,7T2} = (2c(M,iTi, 7 ^ 2 )- Let 

Mx{l) , a,(l)- 5 (l) 


< 


M(l), Y 2 {x) = {D^K^){ 1 ,x). 

(31) 


Proof: We start the proof by observing that since 
{M,Ki,K 2 } S as per Theorem [T] these polyno¬ 

mials define a positive operator V — X^m,Ki,K 2 } ^uch that 
e||tt’|p < {WjVw) < for some 0 > 0. Furthermore, 

by Theorem [2 there exist bounded and continuously differen¬ 
tiable functions M, and K 2 such that P~^ = K K } 
satisfying l/0||r(;|p < {w,'P~^vj) < l/e||t(;|p. We now 
propose the Lyapunov function 

V{w) = (fP~^w,w) = (V~^w,PV~^w) 

= f {V~^w){x)M{x){V~^w){x)dx 
Jo 

+ [ [ {'P~^w){x)Ki{x,^){V~'^w){()didx 

Jo Jo 


{V ^w){x)K 2 {x,^){'P ^w){^)d^dx. 


/ 0 J X 


Let y = V~^w. Note that if w G 77^(0,1), then y = V~^w G 
77^(0,1). Now, since l/0||r(;|p < (ru, V~^w) < l/e|ltt;|p, we 
have that the Lyapunov function is upper and lower bounded. 
Now suppose that 

u{t) =Yi{V~^w){l,t) + f Y 2 {x)(V~^w){x,t)dx 

Jo 

= Riw{l,t) + / R 2 {x)w{x,t)dx. 

Jo 

Since M,K^,K 2 G C°“(0,1) and y 2 is polynomial, we have 
that 7?2 G (7°°(0,1). Therefore, as discussed in SectionHIIl the 
closed loop System (l25T l - (|26]) admits a solution w G 77^(0,1) 
which implies y = 'P~^w G 77^(0,1). Again, the calculation 
of the time derivative V and its reformulation is lengthy. It 
involves integration by parts, the Wirtinger inequality and 
the assumption a{x) > a. This proof is in the appendix as 
Lemma |4] which establishes that for any w which satisfies 
Equations (l25l l - (l26l l. 


2 2 a(l) 

If the control input u{t) is defined as 

u{t) =Yi{V~^w){\,t) + f Y 2 {x){V~^w){x,t)dx 

Jo 

=Riw{l,t) + / R 2 {x)w{x,t)dx, (32) 

Jo 

where P~^ is as defined for V = X[m,Ki,K 2 } Theorem\^ 
and X{m,Ki,K 2 } as defined in Q, then there exists a scalar 
7 > 0 such that for any initial condition w{0) G P (where 
P is as in Equation the solution w of (I25l l - (I261 l exists, 
belongs to c»), [ 0 , 1 ]) and satisfies 


lk(7)ll < 7lk 


t > 0. 


Viw(t)) < {y{t) 

’ ^{M,Ki,K2}y(^')I 

+ y{l,f)Ny{l,t) + 2y{l,f)a{l)M{l)yx{l,f), 

where N = a(l)Ma;(l) -|- (&(1) — aa;(l))M(l). Now, 
since {—M — 2pM,—ki — 2pKi,—k2 — 2 ^X 2 } G 
5 di.d2.o, we have that < -2/rT’ and 

hence (y(7), ^ -2p {y{t),Vy{t)) = 

—2p {w{t), Applying this to the inequality, we get 

V(w{f)) < — 2p {w(t), 

+ y{l,t)Ny{l,t) + 2y(l,t)a{l)M{l)yxil,t). 

A sufficient condition for stability, then, is that 
2 ?/(l)a(l)M(l)?/2^(l) < —y{l)Ny{l). Unfortunately, 

however, our control input enters via uix(l) and not j/x(l)- 
To see the relationship between Wx{f) and yx{L), we expand 
the former and then solve for the latter as follows 

Wxil,t) =Mx(,l)yil,t) + M{l)yxil,t) 

+ [ (DiXi)(l,x)y{xR)dx, (33) 
Jo 
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where solving for M{l)yx{l) yields 

=Wx{l,t) - M^{l)y{l,t) 

- [ {DiKi){l,x)y{x,t)dx. ( 34 ) 
Jo 

This implies that the Lyapunov function satisfies 

< - 2yV{w{t)) +y{l,t)Ny{l,t) 

+ 2y(l, f)a(l) (wx(l, t) - M^{l)y{l, t)) 

-2y{l,t)a{l) f {DiKi){l,x)y{x,t)dx. 

Jo 

Now, examining the proposed controller, we obtain 

= u{t) = Riw{l^t )/ R2{x)w{xR)dx 

Jo 

= Ri{Vy){lR)-\- [ R 2 {x){Vy){x,t)dx 

Jo 

= Fiy(l,t)+ [ Y2{x)y{xR)dx, 

Jo 

which is expressed in the new optimization variables Yi and 
Y 2 . Now, plugging ■u;a;(l) = Yiy{l) + fg Y 2 (x)y(x)dx into 
the time-derivative of the Lyapunov function, we get 

V(w(t)) < — 2fiV(w(t)) 

+ y(l, tf(N + 2a(l)Yi - 2a(l)M,(l)) 

+ 2y{l,t)a{l) [ Y 2 {x)y{x,t)dx 

Jo 

-2y{l,t)a{l) f {DiKi){l,x)y{x,t)dx. ( 35 ) 
Jo 

By inspection, we see that the stability conditions are now 

N + 2a{l)Yi - 2o(l)M^(l) < 0 

and Y 2 (x) = x). This then implies that V{w{t)) < 

—2y,V{w(t)) for all f > 0 and hence V{w{t)) < 

T^(w(0))e“^^‘. Since < 9V{w), we have 

||w(f)|| < \j9V(w{t)) < yj9V{w(Qi))e~^* 

< vW^||u;(0)||e-^‘. 


At this point, it is significant to note that given values for the 
variables Yi, Y 2 , M, Ki and K 2 , the controller gains i?i and 
i ?2 can be found by calculating M, and K 2 via Theorem 
[3 and using the formula 

+ / R2{x)w{x)dx 

Jo 

= Yiy{l)+ f Y 2 {x)y{x)dx 

Jo 

= Yi{V~^w){\) f Y 2 {x){V~^w){x)dx 
Jo 

= yiM(l)u>(l) + j YiYi{l,x)w{x)dx 
Jo 

+J y2{x)(^J Yi{x,^)w{^)d^dx+J K^ix, ^)w{^)ddj dx 


= YiM(l)r(;(l) + f YiYi{l,x)w{x)dx 
Jo 


Y 2 i 0 yii^,x)d^+ / l" 2 ( 0 K 2 (^,a;)d^ w{x)dx 


^0 \^x Jo 

so that 

i?i =yiM(l), (36) 

R2ix)=YiK,{l,x)+ [ Y2{mii^,x)d^+f Y2{0K2i^,x)d^, 

Jx Jo 

(37) 

where we have used the identity 

f{x,^)didx= [ [ f{x,i)dxd^= [ [ f{^,x)d^dx, 


0 Jx 


0 ^0 


0 ^0 


and the fact that =K 2 (C,a;). 


IX. Observer Synthesis 

In Section IVIIII we described LMI conditions under which 
one can obtain controller gains Ri and R 2 {x) such that input 
u{t) = Riw(l,t) + fg R 2 (x)w(x,t)dx ensures exponentially 
stability of the closed-loop system. However, this form of 
controller requires measurement of the state w{x,t) at every 
point X G [0,1] at all times. Implementation of such a 
controller is problematic as such a distributed measurement 
is unlikely to be readily available. A more common scenario 
is one in which we may only measure the value of w{x,t) 
at discrete points in the domain. In particular, we assume 
that only a single measurement is available at the boundary 
of the domain, and in particular, at v{t) = 111 ( 1 ,t). As 
discussed in Section mil in this scenario, we seek to find 
an estimator/observer which will yield a real-time estimate 
of the state of the system at all points and which, if used in 
closed-loop, will ensure exponential stability of the closed- 
loop. Specifically, our observer is a dynamic system with 
input v{t) = w(l,t) and output w(x,t), where w{x,t) is 
the estimate of the state at time t. We adopt the Luenberger 
observer framework discussed previously, which implies that 
the dynamics of the observer are given by 

lit (x,t) = a(x)wxx (xyt) + b{x)wx (x, t) 

+ c{x)w{x, t) + Li(x) (v(t) — v{t)), (38) 

w{0,t)=0, Wxil,t) = u{t) + L 2 (v(t) - v{t)), (39) 


where v(t) = w(l,t) is the predicted output and the scalar L 2 
and function Li(x) are gains which map error in this predicted 
output to the dynamics of the observer state. In the following 
theorem, we seek conditions on Li and L 2 which ensure that 
if i?i and i ?2 are as defined in Theorem |4] and the controller 
is defined as 


u(t) = Riw{l,t) -f 


R 2 (x)w(x, t)dx, 


(40) 


Jo 

then Equations (l3Ft - ( |39] ) coupled with Equations ( l25l l - ( l26b 
and Equation ( l40l ) define an exponentially stable system. 

Our approach is based on the separation principle 13 
Chapter 5], ITSl Chapter 5]. Specifically, we consider the error 
dynamics of the PDE coupled with the observer dynamics 
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in Equations ( l38l l - (l39l l. That is, if we define the error as 
e = w — w, then this quantity satisfies 

etix, t) = a{x)exxix, t) + b{x)ex{x, t) + c{x)e{x, t) 

+ zi{x,t), (41) 

e(0,f)=0, 63 ,( 1 , f) = Z 2 (f), (42) 

where the feedback signals zi and Z 2 are defined as 

zi{x,t) := Li{x)e{l,t) and Z 2 it) := L 2 e{l,t). (43) 

The key point is that the error dynamics do not depend on 
the choice of controller gains Ri and i? 2 - In the following 
theorem, this will allow us to choose observer gains Li and 
L 2 which stabilize the error dynamics. Then, in Theorem | 6 ] 
we will show that if the controller gains are chosen as 
per Theorem |4] and the observer gains are chosen as per 
Theorem |5] then the coupled dynamics are stable in both the 
state and state estimate. Unlike for controller synthesis, the 
conditions for stabilization of the error dynamics are based 
on the use of a simple Lyapunov function V (e) = (e, Ve) 
where the operator V = defined by some 

{M,Ki,K2} G Sdi.da.e- 

The following theorem is motivated by the LMI approach 
as defined in Section |IV] wherein as before the variables M, 
Ki and K 2 are defined by a positive definite matrix P and the 
observer variables are scalar T 2 and polynomials Ti and T 3 
(defined by their vector of coefficients). Referring to the LMI 
motivation, these observer variables are similar to the matrix 
T and the observer gains are then recovered as L 2 = P~^T 2 
and Li = +T3), which is similar to the LMI observer 

gain matrix L = P~^T. 

Theorem 5. Suppose there exist scalars e,S > 0, 

di, d 2 , di, d 2 G N and polynomials M, Ki and K 2 such that 

{M,Ki,K 2 } GEdi,d 2 ,e, and K 2 i 0 ,x)= 0 . 

Further suppose 

{-M - 26M, -k, - 2dKi, -k 2 - 2dK2} G ^dM^ 

where {M,ki,k 2 } := Ki, K 2 ). Let ^ Ki and K 2 

define K K } Theorem\^and 

L 2 :=(a(l)M(l))-iT 2 , (44) 

Li{x) ■=M{x){Ti{x) +T 3 {x)) 

+ r Ki{^k){TiiO+T3md^ 

Jo 

+ [ K^{x,0{Tik) +T3{0)df,. (45) 

J X 

where 

Ti{x) = -0.5((&(1) - akl))Ki{l,x) - a(l)(Z?iiTi)(l,x)), 

(46) 

T 2 < -0.5((&(1) - a,(l))M(l) - a{l)M,{l)) (47) 

T 3 {x) = -L 2 a{l)K 3 {l,x), (48) 

and X{m,Ki,K 2 } and ^ are as defined in Q. Then 

for any e which satisfies ' (14It - (HU with initial condition 


e(0) G T>e (See Equation (I121 l ). there exists a scalar 7 > 0 
such that 

||e(f)||< 7 l|e( 0 )||e-^ f > 0 . 

Proof: We start by observing that since {M,Ki,K 2 } G 
per Theorem[T] these polynomials define a positive 
operator V = X[m,Ki,K 2 } such that e||r(;|p < {w,Pw) < 
0||w|p for some 6* > 0. Lurthermore, by Theorem[3 there exist 
bounded and continuously differentiable functions M, and 
K 2 which define the positive operator K K }■ 

Therefore, since Li G C'°“(0,1), we have that the closed-loop 
error dynamics (HTt - (1421) admit a local in time solution e for 
any Cq G T>e. 

We now propose the Lyapunov function 

Vie)={e,Pe) 

e{x)M{x)e(x)dx 

+ f f e(x)Ki{x,f)ek)dfdx 

Jo Jo 

+ [ [ e{x)K2{x,^)e{f)d^dx. 

Jo Jx 

The derivative of this Lyapunov function is identical to the 
one in Theorem |2] except for the presence of the terms Zi and 
Z 2 defined in (|4^ . Specifically, we have 

14(e) < ^e, 

+ 2 {Ve,zi)p 2 z 2 a{l){Ve)(l) 

+ 2(e,7Wr?3e(l))+e(l)i?2e(l), 

where Ri{x) = (6(1) - axk))Ki{l,x) - a(l)(DiKi)(l,x) 
and i ?2 = (6(1) — aa,(l))M(l) — a(l)M 3 ,(l). In the proof of 
Theorem 121 we had zi = 0 and Z 2 = 0 and hence the stability 
condition was that Ri = R 2 = 0 and that < 

—26P. Lor the observer, we similarly require X^j^ < 

—25V. However, we now have the observer gains Zi{x) = 
Li(a;)e(l) and 63 ,( 1 ) = Z 2 = 7 ^ 26 ( 1 ) which the algorithm 
can choose in order to cancel out i?i and i? 2 - Unfortunately, 
however, these gains depend on M and Ki and the gains 
are currently bilinear with the operator variable V (and the 
functions M Ki, and K 2 which define it). Hence we would 
like to perform a variable substitution. This is complicated, 
however, by the fact that there are two observer gains - one at 
the boundary and one directly injected into the dynamics. Let 
us first examine the second gain Z 2 = U 2 e(l) which appears 
in the term 



Z2a{l){Ve){l) = e{l)L2a{l){Ve){l) 

= e(l)L2a(l)M(l)6(l) 

^ ^ > 

T 2 

+ / e{l)L2a{l)Ki(l,x)e(x)dx 

Jo 

= 6(1)T26(1) + f e{l)L2a{l)Ki{l,x)e{x)dx, 

Jo 

where we have made the variable substitution T 2 = 
L 2 a(l)M(l) which implies T 2 is a scalar variable. The 
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variable L 2 is thereby partially eliminated from the search. 
However, since a(x) >0 and M(x) > 0, given T 2 , the gain L 2 
can later be recovered as L 2 = (a(l)M(l))“^T 2 . Of course, 
this variable substitution has not completely eliminated the 
original variable ^ 2 - To completely eliminate L 2 will require 
assistance from the second gain Li. To see how this is done, 
we examine the second term in which zi appears 

{'Pe,zi) = le, e(l)\=(e, 7 WTie(l))+(e, 7 WT 3 e(l)). 

J^Ti + A^T 3 

Here we have defined a new variable Ti{x) which is de¬ 
fined by Ti(x) := M{x)Li{x) + Ki{x,^)Li{^)d^ + 

K 2 ix,^)Li{^)d^ - Tsix) for which MtiC = TMl^c- 
J^T^c for any c S K. where T 3 will be defined shortly. 
Furthermore, for any T 3 , the map Li ^ Ti is invertible with 

Li{x) :=M{x)(Ti{x) + T 3 {x)) 

+ [\ii^,0{Ti{0+T3{0)d^ 

Jo 

+ [ K2(a:,C)m(e)+T3(e))de 

J X 

if 'P~^ = <T^]y[ K K eliminate the variable 

Li and replace it with Ti and T 3 . The next step, then, is to 
choose T 3 so as to cancel the remaining term which contains 
L 2 . This is done using (e, which we expand to get 

{e,MT^e(l))= [ e{x)Ts{x)e(l)dx, 

Jo 

which we would like to use to eliminate 
/g e{l)L 2 a{l)Ki{l,x)e{x)dx. Clearly, then, the appropriate 
choice for T 3 is 

Tsix) = -L 2 ail)Kiil,x). 

Note that the dependence of T 3 on L 2 is admissible because 
T 3 is not a free variable and L 2 is computed directly from T 2 . 
This means that once feasible values for Ti and T 2 have been 
found, we then calculate L 2 from T 2 , then use L 2 to calculate 
T 3 and then use Ti and T 3 to calculate the gain Li. 

Concluding the proof, the time-derivative of the Lyapunov 
function becomes 

V(e) < (^e, 

+ 2 (e, {A4ti + Mri) e(l)) + e(l) ( 2 T 2 -|- R 2 ) e(l). 

Therefore, if Ti = —Ri, 2 T 2 -|- i ?2 < 0 and {—M — 
26M, —Ki — 2 dKi,—K 2 — 26 K 2 } G p, we have that 


Theorem 6. Given di,d 2 ,di,d 2 G N and 0, suppose 

that there exist polynomials M, N, Ki, K 2 , Si and S 2 such 
that 

{M, Ki,K2} G 

(49) 

{-M - 2pM,-ki - 2^iKl,-k2 - 2pK2} G g, 

(50) 

iT2(0,a;) = 0, (51) 

{N, Si, 82 } G Edi^d 2 ,e, 

(52) 

{_7V - 25N,-k - 25Si, -S 2 - 2,552} G g, 

’(53) 

52(0,a;) = 0, (54) 

where {M,ki,k 2 } = Ki, K 2 ), {N, 81 , 82 } = 

ris^N, Si, 82 ) and 2di and 2 d 2 + 1 are the degrees of M, 
N and Ki, Si, respectively. 

Then, there exist gains Ri, R 2 {x), Li{x) and L 2 such that 

if ^ 

u{f) = Riw{l,t) + / R 2 {x)w{x,t)dx, (55) 

Jo 

and w satisfies Equations (I251 l - (1261) and w satisfies Equa¬ 
tions (1381 ) - (1391 ) with a zero initial condition then 

for some 7 > 0 and any k satisfying 0 < k < min{^, ^}. 

Proof: If the conditions in (09]) - dSB are satisfied, 
then the polynomials M, Ki and K 2 satisfy the constraints 
of Theorem |4] Therefore, we may construct Ri and R 2 {x) 
using (O - (ITtI i. Similarly, if N, Si and S 2 satisfy (l52T i - d^ . 
then the conditions of Theorem |5] are satisfied with M = N, 
Ki = 5i and K 2 = 82 - Thus, we can construct observer gains 
Li{x) and L 2 using (04]) - (05]). Now, let Vc = J^{m,Ki,K 2 }^ 
kc = “ ^{N,Si,S 2} Po = 

Therefore, the theorem conditions imply that Pc < —2p,'Pc 
and Po < —26Vo. 

Using the proof of Theorem |5] there exists a scalar /3i >0 
such that 

Vo{e) <-25Vo{e) - Pie{lf, (56) 

where Vo{e) = {e,Voe). Similarly, for the observer dynamics 
in ([38]) - ( [39] ) with the input (l55l ). using the proof of Theorem0] 
one can prove that there exists a scalar /32 > 0 such that 


U(e) < -2SV{e), 

which, in a similar manner as Theorem|2]establishes exponen¬ 
tial stability of the error dynamics with decay rate S. ■ 

X. An LMI condition for Output-Feedback 
Stabilization 

In this section we briefly summarize the results of the 
paper by giving an LMI formulation of the output-feedback 
controller synthesis problem. 


Vciw) < - 2p.Vc{w) 2 (y, Lie(l)) -b y(l)( 2 a(l)L 2 )e(l) 
-^ 2 y(l)^ (57) 

where y = Vf^w and Vc{w) = (w,Vf^w) = {y,Vcy). 
From (l56] ) - (fSTt we infer that for any r > 0 we have 


rVo{e) + Vc{w) < -2r5Vo{e) + 


y 


y 

\ 

y{f) 

,u 

y(i) 

) 



_e(l)_ 

/ 


( 58 ) 
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where 


U = 


-2fj.Vc 

★ 

'k 


0 Li 

—(32 a(l)L2 , 
* -r/3i 


and Mupad. To help with understanding this process, we define 
several subroutines which perform specific relevant tasks and 
combine them in the pseudo code which would be used to 
obtain the observer-based controllers. 


and the inner product is defined on ^2(0,1) x R x R. Now, 
for any 0 < k < min{i5,/r}, if we choose r > 0 sufficiently 
large, it follows that U < diag(—0,0). Thus, from (l59l l 
we get that 

rVo(e) + Vc(w) < — 2rSVo(e) — 2kVc(w) 

< -2k (rVo(e) + Vc(w)). 

Therefore defining V(w, e) = rVo(e) + Vc{w), we get that 
V(w,e) < —2KV(w,e). 

Integrating in time, 

r (e, Voe) + (w, (ti;(0), VowiO)), (59) 

where we have used the fact that zD(0) = 0 and thus e(0) = 
—tu(0). Now, as discussed, there exist scalars 61,62 >0 such 
that 

e||ef < {e,roe) < 6i\\ef, ^||wf < {w,r-^w) < i||w|| 

t'2 e 

Therefore, using (l59l l we get 

where a = min(re, 1/02)- Thus, 

/ vOi 

||e||,|N|<y^|k(0)||e-«‘. 

Finally, using the fact that ||ir;|| < ||w|| + ||e|| produces 


[M, Ki, K2 ] =mult_semisep (e) 

• Declares polynomial variables M, Ki and K 2 and 
enforces the constraint {M, Ki,K 2 } G 'E.d,i,d 2 ,e- 

[M, Ki, K2 ] =omega_primal (M, Ki, K2 ) 

• Constructs M, Ki and as defined by the map fig 
in (l20l) . 

[M, Ki, K2 ] =omega_dual (M, Ki, K2 ) 

• Constructs M, Ki and K 2 as defined by the map flc 
in (|22l). 

eq_constr(F) 

• Given a set F of univariate/bivariate polynomials, 
declares element wise equality constraint F = 0. 

[M, Ki, K2 ] =inv_op (M, Ki, K2 ) 

• Given {M, Ki, K 2 } = 'E.d^^d 2 ,e calculates the inverse 
multiplier M and kernels and K2 by approximating 
U{x) by performing the integration in (l24l) a finite 
number of times and using a Chebyshev series approx¬ 
imation of M{x)~^. 

[Ri, R2 ] =controller_gains (M, Ki, Ki, M, , K2 ) 

• The function defines Yi and Y 2 {x) using (EB. Conse¬ 
quently, i?i and i?2 {x) are defined using (l36T l and dTTl) . 
respectively. 

[Li, L2 ] =observer_gains (M, Ki, Ki, M, Ki, K2 ) 

• The function constructs T 2 using 63 and sets L 2 us¬ 
ing (I44I 1. Then the function constructs Ti{x) and T^{x) 
using (|46 T i and (l48T l and constructs Li(a;) using (EB. 

A pseudo code for the SOSTOOLS implementation of the SDP 

is presented in Algorithm [T] 


The variables in Theorem |6] are polynomials which are 
parameterized by vectors of coefficients associated to a prede¬ 
termined monomial basis. There are two types of constraints 
on these variables; equality constraints between polynomials; 
and constraints of the form S To test the condi¬ 

tions of Theorem |6l these variables and constraints must be 
converted to a form recognized by an SDP solver such as 
SeDuMi llJTlI . Many of these tasks have already been auto¬ 
mated in SOSTOOLS 1291 and our extended toolbox, Delay- 
TOOLS ll28l . Specifically, SOSTOOLS has functionality for 
declaring polynomial variables and enforcing scalar equality 
constraints. Furthermore, DelayTOOLS ll28l allows the user 
to declare matrix-valued equality constraints and create new 
polynomial variables which satisfy G Sdi,d2,£- Furthermore, 
the multipoly toolbox allows one to manipulate polynomial 
variables in order to construct new dependent polynomials 
such as in {M, ATi, K2} = Ki, K 2 ). Once all variables 

and constraints have been declared, SOSTOOLS converts all 
constraints and variables to a format which can be accepted 
by SDP solvers such as SeDuMi, SDPT3 or MOSEK. The a- 
posteriori polynomial manipulations such as operator inversion 
can be performed using a combination of the multipoly toolbox 


XL Numerical Results 

In this section we test the conditions of Theorems |3 El and |5] 
by applying them to two parameterized instances of scalar 
parabolic PDEs. The first instance is a variation of the classical 
isotropic heat equation. Because this system is well-studied, 
we are able to compare our results with a number of existing 
results in the literature. The second system is an anisotropic 
PDE with arbitrarily chosen coefficients. Both instances have 
an instability term, parameterized by an instability factor, 
A. Eor both systems, we test stability, find controllers and 
construct observer-based controllers. 

Example 1: Our first system is defined as follows. 

wt{x,t) = Wxx{x,t) + \w{x,t), A G R, (60) 

with boundary conditions 

w{0,t)=0, Wx{l,t) = u{t). 

The output of the PDE is v{t) = Eor u(t) = 0, the 

analytical solution of this PDE is given by 

00 

w{x, t) = '^ {wq, (j)n) ((n^x), 

n—1 
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Algorithm 1: Output-feedback controller synthesis. 

Declaring optimization variables: 

1) [M, Ki, K2 ] =mult_semisep (e) 

2 ) [N, Si, S2 ] =niult_semisep (e) 

Constructing polynomials: 

1 ) [M, Ki, K2 ] =omega_dual (M, Ki, K2 ) 

2 ) [N, Si, §2 ] =omega_primal (N, Si, S2 ) 

Declaring affine constraints: 

1 ) eq_constr ( (-M- 2 pM, -Ki- 2 ^Ki, -K2-2/rK2 ) 
-mult_semisep(0)) 

2 ) eq_constr ( (-N- 25 N, -Si- 2 ( 5 Si, -S2-2i5S2) 
-mult_semisep(0)) 

3 ) eq_constr (K2 ( 0 , x) ) 

4 ) eq_constr (S2 ( 0 , x) ) 

if SOS problem is feasible then 
I Return outputs: M, Ki, K2, N, Si, S2. 

Calculating control gains: 

1) [M,Ki,K2]=inv_op(M,Ki,K2) 

2 ) 

[Ri, R2 ] =controller_gains (M, Ki, Ki, M, , K2 

Calculating observer gains: 

1 ) [N, 82,32] =inv_op (N, Si, S2) 

2 ) [Li, L2 ] =observer gains (N, Si, Si, N, ^1, ^2 ) 


d = 3 4 5 6 7 analytic 

A = 0.59 TI9 2:457 2:46 Ziei 21467 

TABLE I; Max. A as a function of di = ^2 = d for which the 
exp. stability conditions of Theorem |2] are feasible, implying 
stability of PDE dbOl l with u{t) = 0. 


where An = A — (2n — l)^7r^/4 and = •\/2sin((2n — 
l)7ra;/2). This implies that Equation dbOl l is unstable for A > 
7r^/4 « 2.467. To test the numerical accuracy of the stability 
conditions in Theorem|2] we found the largest A > 0 for which 
the conditions of Theorem |2] are feasible as a function of the 
parameters di and d2 which dehne the degree of the variables 
M, Ki and 7^2- Table ^presents these results for e,S = 0.001. 
Eor di = d 2 = 7, we can construct a Lyapunov function which 
proves stability for A = 2.461, which is 99.74% of the stability 
margin ^ Ri 2.4674. 

To test the accuracy of the conditions in Theorem |4l we 
hnd the largest A for which the conditions of Theorem |4] are 
feasible with e = 0.001 and fi = 0.001, thereby implying 
the existence of an exponentially stabilizing state-feedback 
controller. Table |II] presents this maximum A as a function 
of the degree di = c?2 = d. The results suggest that for 
sufficiently high degree, a static state-feedback controller can 
be constructed for any value of A > 0. 

To test the accuracy of the conditions of Theorem |5] we 
find the largest A for which the conditions of Theorem |5] are 
feasible with e = 0.001 and 6 — 0.001, thereby implying 
the existence of an exponentially stabilizing dynamic output- 
feedback controller with output v{t) — Table UII] 

presents this maximum A as a function of the degree di = 
d2 = d. The results suggest that for sufficiently high degree. 


d = 7 8 9 10 11 

A = 14.3982 17.9626 22.8645 23.3093 27.1179 

TABLE II: Max. A as a function of di = d 2 = d for which 
the conditions of Theorem |4] are feasible, thereby implying 
the existence of an exp. stabilizing state-feedback controller 
for PDE dlOll. 


d = 7 8 9 10 11 

A = 14.5233 17.7643 23.4406 24.7772 27.8820 

TABLE III: Max. A as a function of di = d2 = d for which 
the conditions of Theorem|5]are feasible, thereby implying the 
existence of an exp. stabilizing output-feedback controller for 
PDE (EOjl. 


a dynamic output feedback controller can be constructed for 
any value of A > 0. 

Example 2: To illustrate the versatility of the proposed 
piethod, we next consider the following arbitrarily chosen 
anisotropic system 

Wt{x, t) = a{x)wxx{x, t) -f b(x)wx(x, t) -\- c{x)w{x, t), (61) 

where a{x) = x^ — x"^ + 2, b{x) = 3x'^ — 2x and c{x) = 
—-f 1.3a:^ — 1.5a; -f 0.7 -f A with A G K.. Although the 
analytical solution to this PDE is not readily available, we 
may use a hnite-difference scheme to numerically simulate 
the system and thereby estimate the range of A for which 
the PDE (ImT i is stable. Specihcally, we hnd that the system 
is unstable for A > 4.66. To determine the accuracy of the 
conditions of Theorem |2] we hnd the largest A for which the 
conditions of Theorem |2] are feasible. Table IIVI lists the largest 
such A using e,S = 0.001 as a function of polynomial degree 
di = d 2 — d. The maximum A for which we can prove the 
exponential stability for is A = 4.62, which is 99.14% of 
the predicted stability margin of 4.66. The < 1% discrepancy 
may be due to conservatism or inaccuracy in the predicted 
maximum A on account of inaccuracy in the discretization or 
poor choice of initial conditions in the simulation. 

To test the accuracy of the conditions in Theorem |4] we 
again hnd the largest A for which the conditions of Theorem |4] 
are feasible with e = 0.001 and p, = 0.001, thereby implying 
the existence of an exponentially stabilizing state-feedback 
controller. Table |V] presents this maximum A as a function 
of the degree di — d 2 = d. The results suggest that for 
sufficiently high degree, a static state-feedback controller can 
be constructed for any value of A > 0. 

To test the accuracy of the conditions of Theorem |5l we 
again hnd the largest A for which the conditions of Theorem |5] 
are feasible with e = 0.001 and 5 = 0.001, thereby implying 


d = 3 4 5 6 7 simulation 

A = 4.37 461 461 462 462 466 

TABLE IV: Max. A as a function of di = ^2 = <^ for which the 
exp. stability conditions of Theorem |2] are feasible, implying 
stability of PDE (ImT i with u{t) = 0. 
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d = 4: 5 6 7 8 

A = 19.0216 36.1359 39.7247 43.5974 44.5219 

TABLE V; Max. A as a function of di = ^2 = d for which 
the conditions of Theorem |4] are feasible, thereby implying 
the existence of an exp. stabilizing state-feedback controller 
for PDE dlB- 


d = 4 5 6 7 8 

A = 18.3090 36.0199 38.0478 40.5931 44.079 

TABLE VI: Max. A as a function of di = d 2 = d for which 
the conditions of Theorem|5]are feasible, thereby implying the 
existence of an exp. stabilizing output-feedback controller for 
PDE dSB- 


the existence of an exponentially stabilizing dynamic output- 
feedback controller with output v{t) = w{l,t). Table |Vl] 
presents this maximum A as a function of the degree di = 
d 2 = d. The results suggest that for sufficiently high degree, 
a dynamic output feedback controller can be constructed for 
any value of A > 0. 

We conclude with the conjecture that the proposed method 
is asymptotically accurate in the sense that, for any A > 0, 
if the PDE O - Q is stable in the autonomous sense, then 
the conditions of Theorem |2] will be feasible for sufficiently 
high di and d 2 - Moreover, we conjecture that if the system 
is observable and controllable for some suitable definition 
of controllability and observability, then the conditions of 
Theorems |4] and |5] will be feasible for sufficiently high di and 
d 2 - We emphasize, however, that this is only a conjecture and 
additional work must be done in order to make this statement 
rigorous and determine its veracity. A further caveat to these 
results is the observation that the maximum degree di and d 2 
for which the conditions can be tested is a function of the 
memory and processing speed of the computational platform 
on which the experiments are performed. Specifically, the 
number of optimization variables in the underlying SDP prob¬ 
lem is determined by the number of polynomial coefficients 
which scales as 0{d^). To illustrate, all numerical experiments 
presented in this paper were performed on a machine with 
8 gigabytes of random access memory, which limited our 
analysis to a maximum degree of di = d 2 = 11 for PDE (l60l l 
and di = d 2 = S for PDE (IhTl l. 

In the following subsection, we illustrate the controllers and 
observers which result from feasibility of the conditions of 
Theorems |4] and |5] using numerical simulation. 

A. Numerical Implementation of Observer-Based Controllers 

To illustrate the observer-based controllers which result 
from feasibility of the conditions of Theorems |4] and |5] we 
take the anisotropic PDE (IhTT i with A = 35. This value of A 
renders the autonomous system unstable. We then synthesize 
controller and observer gains using the results of Theorems E] 
and |5] for di = d 2 = 6, along with the inverse state 
transformation defined in Theorem |3] Eor the inverse state 
transformation, is approximated using a sixth order 



d 

= 1 

2 

3 

4.. .9 

10 

Ki,K2fO 

Ex. 1 

A = 

: 3.91 

4.78 

4.88 

4.88 

4.88 

27.1179 

Ex. 2 

A = 

: 3.51 

7.03 

8.59 

8.59 

8.59 

44.5219 

TABLE 

VII: 

Re-evaluation of the results of Tables |II] and |V] 

with added constraint Ki 

= K 2 

= 0. 




d 

= 1 

2 

3 

4.. .9 

10 

K-i,K2 A 0 

Ex. 1 

A = 

: 3.89 

4.79 

4.88 

4.88 

4.88 

27.8820 

Ex. 2 

A = 

: 3.51 

7.12 

8.43 

8.43 

8.43 

44.079 


TABLE VIII: Re-evaluation of the results of Tables Hill and I VII 
with added constraint Ki = K 2 = 0. 

Chebyshev series approximation and 5 iterations are used to 
define C/qo — Us- The controllers are then applied to the state 
and estimator dynamics, which are then discretized using a 
trapezoidal approximation. The initial state is set to 

(x-0.3)^ (x-Q.7)^ 

Wo{x) = e 2(0.07)^ _ e 2(0.07)2 ^ 

while the initial observer state is set to w{x, 0) = 0. Eig- 
ures |2(a)| - |2(c)| illustrate the state evolution of the system, 
observer and the control effort respectively. Einally, Eigure |3] 
illustrates the integral control gain R 2 {x). Note that its behav¬ 
ior at the boundaries is logical since at x = 0 , the boundary 
condition t) = 0 ensures that no control effort is required. 
Whereas, at a; = 1, the control exerts maximum effort. 

XII. Necessity of Semi-Separable Kernels in the 

Lyapunov Eunction 

Recall that the Lyapunov functions used in Theorems |2l E] 
and 13 all have the form 

V{w) = / w{x)M{x)w{x)dx 

Jo 

Ki{x,^)w{^)d^ + J K2{x,^)w{^)d^^dx. 

As mentioned previously, this form is atypical in the study 
of parabolic PDEs and the reader may question the neces¬ 
sity of the terms Ki and K 2 as their presence significantly 
complicates the analysis and increases the complexity of the 
stability conditions. Therefore, to illustrate the necessity of 
including these terms, in this section we repeat the numerical 
examples presented previously with the added restriction that 
Ki — K 2 = 0 (which translates to Pij — 0 for z ^ j ^ 1 in 
Theorem[T]|. Table lVIll illustrates these results for the controller 
synthesis conditions of Theorem El using the same method¬ 
ology as described in the previous section. These numerical 
tests indicate that while inclusion of Ki and K 2 allows us 
to control the PDE for any A > 0, when Ki = K 2 = 0, 
our method will fail for some A, regardless of the polynomial 
degree di = di = d. As indicated in Table IVIlH the results are 
similar for the observer synthesis conditions of Theorem |5] 

XIII. Comparison With and Relation to Existing 

Results 

In this section, we compare our numerical results with 
several results in the literature which can be used for stability 
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(a) Evolution of closed-loop state w{x^t). (b) Evolution of closed-loop state estimate w{x,t). 


(c) Control input = u{t). 


Fig. 2: Evolution of closed loop system for Example 2 with A = 35 using controller from Theorem 0] and observer from 
Theorem |5] 


analysis and control, including those based on Sturm-Liouville 
theory and backstepping. 


A. Static Controllers Using Stunn Liouville Theory 

The output feedback controllers we construct are dynamic 
in that they rely on an auxiliary set of estimator dynamics 
which must be simulated in real-time. By contrast, static output 
feedback controllers do not use an estimator and instead rely 
only on a gain of the form, e.g. u{t) = —Kv{t) = —Kw{l,t). 
Unfortunately, even for finite-dimensional systems the problem 
of static output feedback design is unsolved when B ^ I. 
That is, there is no LMI or polynomial-time algorithm which 
is guaranteed to find a stabilizing output feedback controller 
if one exists JSl, 113. However, there are numerous results 
which give sufficient conditions for the existence of such a 
controller, often based on the use of a fixed Lyapunov function. 
Eor the parabolic PDE which we consider, Sturm-Liouville 
theory ifTOl Chapter 2] can be used to express conditions for 
existence of static-output feedback controllers. Specfically, for 
u{t) = the stability of ([T]i - (|2]l depends on the 

first eigenvalue of the following Sturm-Liouville eigenvalue 
problem 

^ \ + q{x)w{x) = ^ia{x)w{x), (62) 

ax \ ax J 


where // is the eigenvalue and 


p{x) 



q{x) = c{x) 


P{x) 

a{x) ’ 


cr(a:) 


P{x) 

a{x) 


The boundary conditions for this eigenvalue problem are 
w(0) — 0 and Wx(l) + kw(1) = 0. Eor our system, using 
the properties of the coefficients a(x), b(x) and c(x) it can 
be established that p is continuously differentiable, q and a 
are continuous and there exist scalars pg and uq such that 
p(x) > Po > 0 and a{x) > o-g > 0. If pi is the first eigen¬ 
value of ( l62b . then it can be established using the Rayleigh 
quotient that where /i“ is the first eigenvalue of 

the following constant coefficient Sturm-Liouville eigenvalue 
problem 


Po 


cPw{x) 

w(x) 


qiw{x) 


pC'^aiw{x), 


(63) 


subject to the boundary conditions w{Q) = 0 and tUa;)!) -f 
Kr(;(l) = 0 and where qi and cti are scalars such that 

q{x) < qi and cr{x) < ai. 

Now let us first consider Numerical Example 1, as defined 
in Equation (l60l l in Section |Xl] In this case, we have that 
Po = 1, qi = \ and cti = 1. Therefore, estimating the first 
eigenvalue of (l63t we get that Ri A—tt^. Since, for stability 
we require < 0, for a large enough k > 0, a control 
input of the form u{t) = —can stabilize dhOl l for 

A < TT^. This result is significantly more conservative than 

the results described in Tables HHHIII which yield a stabilizing 
controller for at least A < 27.1179. Of course this result is not 
particularly surprising, as static output feedback controllers are 
a subset of dynamic output feedback controllers. 

Similarly, for Numerical Example 2 (Equation dMT l) we have 
p{x) = — x^ + 2, q{x) = —0.5a::^ + 1.3a:^ — 1.5a; + 0.7 + A 

and a{x) = 1. Thus pg = 50/27, qi = 0.7 -f A and cti = 1. 

Therefore, estimating the first eigenvalue of (l63t we get that 
R! A — 17.58. As before, we require < 0. Therefore 
for a large enough k > 0, a control input of the form u{t) = 
—Kw{l,t) can stabilize dhOl l for A < 17.58. Whereas, from 
Tables IVTlVIl we see that Theorems 0] and |5] yield a dynamic 
output feedback controller for at least A < 44.079. 


B. The Case When A -f A* < 0 


Eor some values of the coefficients a{x), b{x) and c{x) 
we may have that A + A* < 0 on where the operator 
A is defined in 0 and the set Vg is defined in 0. The 
output feedback stabilization of such systems, i.e. systems with 
< 0 and collocated control/observation, is considered 
in 13 ■ The authors in 0 show that for such systems there 
exists a scalar k > 0 (possibly k = oo) such that the control 
u{t) = —Kv{t) exponentially stabilizes the system. We wish 
to see if our methodology offers a performance gain over the 
controller proposed in 0. If we choose a{x) = 1, b(x) = 0 
and c{x) = 7r^/4, then 


A^ 




(64) 


dx'^ 4 

Applying integration by parts and Lemma [T] it can be estab¬ 
lished that A -f A* < 0 on T>g. If we apply a controller of 
























16 


d = 6 7 8 9 10 11 

<5 = 8.01 12.7 17.21 20.31 22.66 25.78 

TABLE IX; Max. exp. decay rate 5 as a function of polynomial 
degree, di = d 2 = d for Equations ([T]l - (Ell with A as in 
Equation ( l64l l for which we can construct output feedback 
controllers using Theorems 0] and |3 


the form proposed in then u{t) = —Kv{t) = 
for some k > 0. Using the theory in Subsection IXIII-Al it is 
easily established that even for an arbitrarily large k > 0, the 
closed loop system state will decay with a rate close to, but 
less then 37r^/4. Whereas, from Table HXl we observe that for 
di = c?2 = 11 we can construct an output feedback controller 
with a minimum exponential decay rate of 25.78, a significant 
improvement over 37r^/4. 



Eig. 3; Control gain R 2 {x). 


C. Backstepping 

Backstepping is a well-known alternative for the construc¬ 
tion of stabilizing controllers for parabolic PDEs. Specihcally, 
the backstepping approach defines a control law which, when 
coupled with an invertible state transformation, converts the 
controlled parabolic PDE to the form of a desired stable 
PDE (the target system). Although backstepping is not an 
optimization-based method and does not explicitly search for a 
Lyapunov-based stability proof, it turns out that the existence 
of a backstepping controller typically implies the existence of 
a Lyapunov function of the Eorm (fTTI l, defined by a multiplier 
M and semiseparable kernels Ki and 7^2 ■ To demonstrate the 
existence of this Lyapunov function, let us consider the system 
dehned by Example 1, 

wt{x,t) = Wxx{x,t) + ^w{x,t), (65) 

ti;(0,f) = 0, Wx{lA) = u{t), (66) 


then for any solution w of Equations 

z{x,t) ='w{x,t) — / E{x,^)w{^,t)d^j 

Jo 

is a solution of the target system in Equations (IhTl i - dhST l. 
Eurthermore, if the map £ : w z is invertible, then stability 
of the target system implies stability of the original closed- 
loop PDE. Eor the example problem given, this E is obtained 
as a solution of a kernel-PDE and can be found explicitly 
as lITSl 


E{xA) = -Xi- 




0 < c < a; < 1, (69) 


where Ii is the first order modihed Bessel function of the hrst 
kind. Moreover, £ has an inverse of the form 


{£ ^z){x) = z{xA)+ [ E{x,^)z{^,t)d^, 
Jo 


(70) 


where 


F{x ,0 = -K 




0 < C < a; < 1, (71) 


where Ji is the first order Bessel function of the hrst kind. 
Using properties of Bessel functions, it can be shown that both 
kernels E and E are bounded on the domain {{^,x) : 0 < 
^ < a; < 1}. This implies that both £ and £~^ are bounded 
with induced norms which we denote by HfUr; and 

Now, to understand how this backstepping transformation 
implies the existence of a Lyapunov function with semi¬ 
separable kernels, we hrst note that stability of the target 
system in Equations (l67l) - (l68l) is established using the simple 
Lyapunov function 

Vtargetiz) = / z{xf dx = {z,z) , 

for which, using (l67l i - (l68l) . integration by parts and Lemma[T] 
we obtain 




(72) 


for any 2: which satishes ( l67T i - (l68l l. where e = ^- This 
implies 

Vtarget{z{t)) < (^(0)) ^ \\z{x,t)\\ < 6“ 5 ‘ || ^(a;, 0) ||. 


Now, for the original system we dehne the Lyapunov 
function 


Vplant{w) = {£w,£w). (73) 


where A > 0. Now dehne the target system 

zt(x,t) = Zxxix,t), (67) 

z{0,t) = 0, 2 x(l,i)= 0. (68) 

The key backstepping result is that there exists a function E 
such that if 

u{t) = E{l,l)w{i^,t) + f {DiE){l, x)w{x,t)dx, 

Jq 


Now, since for any solution, w{t), of the original system, z = 
£w{t) is a solution of the target system, we have that 

^Vpiant{w{t)) = ^ {£w{-,t)Xw{-,t)) 

= ^ {z{-,t),z{-,t)) 

— ^ €.\^target{_z(t)) 

= -eVtarget{^w{t)) = -€ {8w{t), £w{t)) 
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= -eVplant{w{t)). 


Therefore, 


Vplant{w{t)) < e "^VplantiwiO)), 


which means 


\\£w{-,t)\\<e 2*||f.u;(.,0)||. (74) 

Boundedness of 8 . and £~^ now implies IkWII < 
||^"lc||fu'(f)|| and ||fu;(0)|| < ||f||£||'u;(0)||, which yields 

\\wm < \\£-^cm\ce-nwm\, 

which proves that Vpiantiw) = ||fw|p establishes exponential 
stability of the original system. 

We now show that Vpiant{w) has a form consistent with 
Theorem 0] Expanding 


Vplant{w) = {£w,£w) , 


we get 


Vplant{w) = / w{xY - / / w{x)E{x,^)w{^)d^dx 

Jo Jo Jo 

n w{x)E(^, x)w(^)d^dx 

+ / / / w{x)E{^,x)E{^,rj)w{ri)dr]d^dx. 

Jo Jx Jo 

Changing the order of integration twice in the last integral and 
collecting like terms, we obtain 

Vplant{w) = / W{xfdx 

Jo 

+ [ f w{x)Hi{x,^)w{^)d^dx 

Jo Jo 


+ f ( w{x)H2{x,^)w{^)d^dx, 
Jo Jx 


where 

Hiix,^)=f E{f],x)E{r],^)dr] - E{x,0, 

J X 

H2{x,C}=j^ E{T],x)E{r],^)dr] - E{^,x), 

which has the form of a Lyapunov function consistent with 
Equation (fTTI i using a semi-separable kernel where we have 
M{x) = 1, Ki = Ell and K 2 = 772- In a similar manner, if 
we define V = T’{/,Gi,G 2} where 

Gi{x,^)=f E{r],x)F{T],^)dr] + E{x,^), 

J X 

G2{x,i)=j^ E{r],x)F{T],^)dT] + E{^,x), 

then 'P~^ = X{i.Hi,H 2 } and hence 

Vplantiw) = {V-^w, VV-^w ), (75) 

which is a form consistent with Theorem 0] Thus we conclude 
that for this class of systems, if we assume the function F 


may be approximated by polynomials, then the existence of a 
backstepping controller implies the feasibility of Theorem 0] 
for some degree. Despite this similarity, there are, of course, 
differences between the proposed method and backstepping. 
Specifically, our approach is optimization based, whereas the 
search for the backstepping transformation is not. Advantages 
of the proposed method include the ability to analyze stability 
of autonomous PDEs and simple extensions to robust control 
of PDEs with parametric uncertainty via Positivstellensatz 
results llini. 


D. Finite-Dimensional Approximations 

In this subsection we consider the merits of the SOS 
approach with respect to finite-dimensional approximation. 
That is, we consider whether there are advantages over model 
reduction techniques wherein the PDE is reduced to a set of 
coupled ODEs - as in, e.g. m. 

Before continuing, we note that establishing a suitable 
metric for comparison of finite-dimensional and infinite¬ 
dimensional approaches is complicated by the fact that that 
the methods proposed in this paper are suboptimal. That is, we 
are not seeking observer-based controllers which are optimal 
in any sense. Rather, we simply seek observer-based con¬ 
trollers which establish closed-loop stability. In this sense, our 
methods are roughly equivalent to existing finite-dimensional 
approaches in that for all numerical examples considered, we 
are able to construct observer-based controllers for suitably 
high polynomial degree. In a sense, then, one could argue 
that finite-dimensional approaches are superior in that they are 
able to go beyond stabilization and construct optimal observer- 
based controllers using a suitably high level of discretization. 
In practice, however, our experience has shown that there are 
disadvantages to discretization-based methods such as pole- 
placement. Specifically, we have seen that if the reduction 
scheme is not carefully chosen, discretization may result in 
loss of controllability or poorly conditioned controllability 
matrices. To illustrate, consider the following model; 

Wtix,t) = Wa:xix,t) -f 

w{ 0 ,t) = 0 , = u(t). 


One approach to reduction of this PDE to a system of ODEs 
is to use a finite difference method to approximate the spatial 
derivative as 


c (Xj t) 


)(x -f Ax 2 ,t) — w{x, t) 

Aa;i -f Ax 2 V 

2 f w{x,t) — w{x — Axi,t) 


Axi -f Ax 2 


Axi 


where Axi is the step size to the left of x and Ax 2 is the 
step size to the right. Using this scheme we obtain an ODE 
model of the form 


^m(^) ^ (76) 

where and A™ G and m G N is 

the order of reduction. While relatively straightforward, this 
approach creates significant technical challenges. For example: 
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m 5 10 20 

cond{C{A'^, B’^)) Ri lO’ 10*’^ 

TABLE X: Condition number of C(A™, B^) as a function of 
order of reduction m. 


a) Controllability of the Reduced Model: The reduced- 
order model must be chosen so as to maintain the properties 
of controllability and observability. In most cases, however, 
there is no guarantee that a hnite-difference approximation 
scheme will preserve these properties. For example, for the 
hnite difference scheme dehned above, it is known that if 
the original system is controllable and a uniform grid size is 
chosen, then the reduced system is also controllable. However, 
if one were to chose a non-uniform grid, then controllability 
is no longer guaranteed. For example if one were to chose 
a logarithmic grid, for m > 13 the reduced model is not 
controllable (although it is still stabilizable). In such a case, 
the performance of the closed loop system will be limited by 
the location of the uncontrollable eigenvalues. 

b) Ill-conditioned Controllability Matrix: Now suppose we 
wish to perform pole placement by applying Ackermann’s 
formula to the reduced order model. As mentioned, it can be 
shown that the reduced order model in (l76l l is controllable 
for any m G N when derived using uniform step sizes 
(Aa;i = Aa;2) as established by the Hautus test. However, 
the pole placement problem (which is similar to our condition 
for exponential stabilization with desired decay rate) relies on 
inversion of the controllability matrix C{A^, B^) - a step 
which is numerically sensitive to conditioning of C{A^, B^). 
This is problematic since, as seen in Table |2 the controllabil¬ 
ity matrix for this system is ill-conditioned and the condition 
number worsens as the level of disretization m increases. This 
implies that as the level of discretization increases, numerical 
errors may dominate - potentially resulting in unstable or 
unpredictable controllers. Naturally, these issues are well- 
known and have been addressed in the literature through 
methods such as robust place placement ioi or Galerkin 
schemes ifT^ . The advantage of the SOS approach, however, is 
that the controllers are provably stable at the pre-lumping stage 
and thus the only numerical concern is implementation, which 
does not appear to be sensitive to issues such as condition 
number. 

XIV. Alternative Boundary Conditions 

The results of this paper may be readily adapted to other 
types of boundary conditions. Specihcally, the conditions of 
Theorems |2] |4] and |5] can be easily modihed to consider 
alternative boundary conditions. Although economy of space 
prohibits us from presenting these conditions in full, in this 
section we give the results of numerical tests performed 
using Dirichlet, Neumann and Robin boundary conditions. 
Specihcally, for the two PDFs dbOl l and (IbTl l which dehne 
Examples 1 and 2, respectively, in Section IXj we consider 
the boundary conditions and the outputs as listed in Table |Xl] 



Boundary Condition 

Output v(t) 

Dirichlet 

rH(0, t) = 0 
w{l, t) = u{t) 

Wxil, t) 

Neumann 

Wx(0, t) = 0 

UIa;(l, t) = u(t) 

w{l,t) 

Robin 

ui(o, t) -y wxio,t) = 0 
w{l,t) + = u{t) 



TABLE XI; Alternative boundary conditions and outputs for 
PDFs dSnil and (gB- 


d = 8 9 10 11 

Dirichlet A = 17.7634 22.8645 23.3093 27.1179 

Neumann 14.8163 17.1814 21.8781 21.8781 

Robin 13.8367 16.6565 18.6050 18.9758 

TABLE XII: Max. A as a function of polynomial degree, 
di = d 2 = d for PDF dbOl l with boundary conditions as in 
Table IXIl for which we can construct output-feedback boundary 
controllers. 


Tables IXlIl and IXlIIl illustrate the maximum A for which we 
can construct output-feedback based controllers as a function 
of di = d 2 = d for PDFs dbOl l and dSB, respectively, for the 
boundary conditions listed in Table IXIl using exponential decay 
rates of 6 = = 0.001. Similar to the observation made in 

Section IXII the numerical results in this section suggest that 
our methodology is asymptotically accurate for the considered 
alternative boundary conditions, that is, given any A > 0, we 
can construct controllers/observers by choosing a large enough 
di = d 2 = d. A more detailed study of alternative boundary 
conditions can be found in the thesis work of m. 

XV. Conclusion and Future Work 

We have dehned an algorithmic, polynomial-time approach 
to the design of observer-based controllers for a general 
class of scalar parabolic partial differential equations using 
measurements and feedback at the boundary. The results use 
polynomials and semidehnite programming to parameterize 
a convex set of positive Lyapunov functions on the Hilbert 
space L2. By combining these Lyapunov functions with an 
invertible state transformation, we obtain convex conditions 
for stability, controller synthesis and Luenberger observer 
design. Furthermore, we have tested our results using pa¬ 
rameterized numerical examples in order to show that the 
stability conditions are accurate to several signihcant hgures 
and the synthesis conditions yield controllers for a large class 
of controllable and observable systems. Furthermore, we have 



II 

6 

7 

8 

Dirichlet 

A = 36.0199 

38.0478 

40.5930 

44.079 

Neumann 

29.8492 

31.1447 

31.1447 

34.1584 

Robin 

24.6490 

27.8503 

27.8503 

29.4373 

TABLE XIII; 

Max. A as a 

function of polynomial degree, 


di = d 2 = d for PDF (IMT l with boundary conditions as in 
Table IXIl for which we can construct output-feedback boundary 
controllers. 
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adapted the approach to three alternative classes of boundary 
measurements and actuators. Finally, we have performed a 
series of comparisons with existing results in the literature, 
showing, e.g. that the method is analytically equivalent to 
backstepping for controller synthesis and furthermore is nu¬ 
merically competitive for the examples considered. By using 
an optimization-based algorithm defined by polynomials, the 
results presented here have the advantage that they may be 
further extended to the problem of nonlinear stability analysis, 
robust control, and control of coupled, multivariate, hyperbolic 
and elliptic PDFs - topics of ongoing research. 


Appendix 


To facilitate presentation in this appendix, we use the 
following lemmas. The first is simply a restatement of the 
Wirtinger inequality 

Lemma 1 (1311). Let z G 1) be a scalar function. Then 


/•i 4 /"i 

/ (z(x) — z( 0 ))'^dx < ^ / zJxYdx. 
7o Jo 


The second lemma is accomplished by splitting the integral 
in two parts and applying a change in the variable of integra¬ 
tion to the second part. 


Lemma 2. For any bivariate polynomials K and P the 
following identity holds for any w G ^2(0,1) 


K{x,^)w{^)d^ + J P{x,^)w{^)d^'^ dx 
Jo Jo + 

+ /y + K{^,x)]w{^)d^dx. 


Lemma 3 (Analysis). Given polynomials a, b and c with 
a(x) > a > 0, for all x G [0,1], suppose that there exists 
a scalar e > 0 and polynomials M, Ki and K 2 such that 


{M, Ki,K2} G ^di,d2,ei 

(6(1) - a,(l))Ki(l,x) - a(l)(niKi)(l,x) = 0, 

(6(l)-a,(l))M(l)-a(l)M,(l)<0, 

K2(0 ,x)=0 . 


Let 


is self-adjoint, we can write V{w{t)) = 2 {wt,Pw). Using 
Equation (fTsT i we expand this out to get 


V{w{t)) = 2{wt,Pw) =2Y,^„ 


(77) 


where 
Ti = 


/ Wxx{x,t)a{x)M{x)w{x,t)dx, 

Jo 

T 2 = Wx{x,t)b{x)M{x)w{x,t)dx, 

Jo 

r3=V' / Wxx{x,t)a{x)Ki{x,^)w{^,t)d^dx, 
i=i J^^ 

r4=^ / Wx{x,f)b{x)Ki{x,^)w{^,f)d^dx 

i=i J^i 

Ts = / w{x,t)'^M{x)c{x)dx 

Jo 

+ y / w{x,t)c{x)Ki{x,^)w{^,t)d^dx, 

i=i 

where Ai = {(^,ai) : 0 < ^ < x < 1} and A2 = {(C) ■ 

0<x<C<l}- Applying integration by parts twice and 
using the boundary condition w{0,t) = Wx{l,t) = 0 yields 

ri=— f Wx{x,t)'^a{x)M{x)dx 

Jo 


1 T — 

2 ./n 9x2 


[a{x)M (x)] w{x, t)‘^dx 


- - {ax{l)M{l)+a{l)Mxil))wil,t)^. 

Since a{x) > a > 0 and {M, Ki, K 2 } G we have 

a{x)M{x) > ae. Thus, by application of Lemma [T] we get 

TT^ 

— / Wx{x,t)^a{x)M{x)dx < —;rci£ / w{x,t)^dx. 

Jo 4 Jq 


'0 

Therefore, we conclude that 

Ti < 


iy w(x,f)2 [a(x)M(x)] - yae^ dx 


- i(a,(l)M(l) + a(l)M,(l))u;(l,f)^ (78) 

Again, applying integration by parts once and using i(;(0, t) = 

0 , 


V{w) = {w,X^m,KuK 2 }w) , 

where X^m,Ki,K 2 } A as defined in Q. Then, for any w which 
satisfies Equations (El - (El, 

where {M, ATi, 7^2} = Ki, K 2 }. 

Proof: Let V = X{m,Kx,K 2 } so that V{w) = {w,Vw). 
If w satisfies (El - (El, then taking the time derivative 
of V{w{t)) and since {M, Ki, K 2 } G Ed^^d 2 ,e implies P 


r2=-iy w{x,tf^ [6(x)M(x)] dx+^b{l)M{l)w{l,t)‘^. 


(79) 


Since (M, ATi, A:2} £ ‘^di,d 2 ,e, we have Ari(x,C) = K 2 {^,x) 
and thus Ki{x,x) = 7^2(x,x). Exploiting this property, 
the constraint 7^2(0, x) = 0, and the boundary conditions 
w{ 0 ,t) = Wx{l,t) = 0, we apply integration by parts twice 
to obtain 


r,3 = 


'j{x,tf 


dx 


[a{x){Ki{x,0-K2{x,^))] 


dx 


J 
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[a{x)K,{x,^)]^ w(i,t)did 

— w{l^t) / axil)Ki{l, x)w{x,t)dx 

Jo 

— w{l,t) / a{l){DiKi){l, x)w{x,t)dx. 

Jo 

Applying Lemma|2]and using Ki{x,^) = K 2 {^,x), we get 

^ 3=1 w{x,t) 

Jo 


^ where y = 7^ and {M^Ki,K 2 \ &^c{M,Ki,K 2 \. 

Proof: Taking the time derivative of V{w{t)) and since 
V~^ is self-adjoint, we obtain 

V{w{t)) = 2 {wt,V~^w) 

= 2 {a{-)wxx + b{-)wa: + c{-)w,'P~^w) 


o2 Q 


■^Hx){Ki{x,^) - K2{x,^))] 


= 2^r„, 


(83) 


J i = Xy 


/ (aa;(l)A:i(l,x)+a(l)(L>iA:i)(l,a;)) w{x,t)dx 
Jo 


n—1 


where y = V and 


^ - / w{x,t) 


1 _ 

T 

a{x)Ki{x,^) 

1 

_1 


a{^)Ki{x,^)_ 


t)d^dx. 

(80) 


ri=y a{x) — {M{x)y{x,t))y{x,t)dx, 
r 2 =y b{x)-^{M{x)y{x,t))y{x,t)dx, 


Applying integration by parts once and following the same 
procedure as for Fa, we get 

) nT 


10 

2 .1 


92 


r==E/. 


r- = -E 


1 


w{x, t) 


dx 

_a_ 

LaiJ 


b{x)Ki{x,C) 


2=1 

2 


+ u>(l,f) f b{l)Ki{l,x)w{x,t)dx. (81) 

Jo 

Finally, employing Lemma |2] produces 

^ 5 = [ w{x,t)'^M{x)c{x)dx 

Jo 

^' 1 f 

+ ^-y^ w{x,t){[c{x)+c{^)]K^{x,^))w{^,t)d^dx. 


wif,t)d^dx ^4-^/ Kx) ^ 


^ K^ix,^)y{^,t)d ^j y{x,t)dx, 
K^{x,C)y{i,t)d^ ) y{x,f)dx, 


(82) 


dx 

2=1 ■ 

Fs = / c{x)M{x)y{x,ty‘dx 

Jq 

i=i d^i 

where Ai = {(^,x) : 0 < ^ < x < 1}, A2 = {(C, x) : 0 < 
X < ^ < 1}, /3i = [0,x] and {32 = [x, 1]. Before proceeding 
we calculate y{0,t). The definition y = V~^w implies 

w(0,t) = M{0)y{0,t) + [ K 2 ( 0 ,x)y{x,t)dx. 


Finally, we combine the terms (El - dSl into the deriva¬ 
tive (iTTl i and use the constraints 


Therefore, since ?ii(0,f) = 0 and 212(0, x) = 0, we get 
2/(0, f) = 0. Now, since M{x)a{x) > ae and y{0,t) = 0, 
applying integration by parts twice and using Lemma [T] 
produces 

to eliminate extraneous terms, thereby completing the proof. 7i<— f (a XX {x)M{x)+a{x)Mx:x{x)-'^a^ y{x,t)'^dx 


(6(1) - a,(l))ili(l,x) - a(l)(fliili)(l,x) = 0, 
(6(1) - a,(l))M(l) - a(l)M,(l) < 0, 


1 


Lemma 4 (Controller Synthesis). Given polynomials a, 6 and 
c with a(x) > a > 0 , for all x S [0,1], suppose that there 
exists a scalar e > 0 and polynomials M, Ki and K 2 such 
that 


'a{l)Mx{l) - ax{l)M(l))y{l,tf 


Let 


{M, ATi, 112 } G il2(0,x)=0. 


V{w) = {w,V ^w) , 


+ a(l)M(l)2/,(l,f)2/(l,f). (84) 
Similarly, applying integration by parts once yields 
1 

^ 2 = 7 : / {b{,x)Mx;{x)-bx{x)M{x))y{x,tf-dx 

+ i6(l)M(l)2/(l,<)2. (85) 

Applying integration by parts twice and Lemma H yields 


where V = ‘^^d X{m,Kx,k- 2 } defined in (H. 

Then, for any w which satisfies 


r., = 


/ U.) 


^[K,{x,0-K2{x,0\ 




V(w{f)) <{y{t),Xf^M f^^ f.^}y(f)^ 

+ [a(l)M,(l) + (6(1) - a,(l))M(l)] y{l,tf 

+ 2a(l)M(l)2/x(l,f)2/(l,f), 




a(e)^ 


Ki{x,0 

Ki{x ,0 


y{x, t)^dx 

y{^,t)d^dx. 

( 86 ) 
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In a similar manner as for Fa, we obtain 

2 


r4=E 


1 


i=l 


'Ai 


yix,t) 






9? 


K^{x,i) 


yii,t)d^dx. 

(87) 


Finally, applying Lemma E] to Fs produces 

Fs = / c{x)M{x)y{x,tYdx 
Jo 

1 ' f 

+ 2 ^J^ yix,t)(c{x) + c{0)Ki{x,^)y{C,t)d^dx. 


Substituting Equations 
proof. 


into 


( 88 ) 

completes the 
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